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1. Introduction 

1.1. Matrix models per se 

The purpose of these lectures is to present basic matrix models as practical combi- 
natorial tools, that turn out to be "exactly solvable" . In short, a matrix model is simply 
a statistical ensemble of matrices with some specific measure, here given as an invariant 
weight, to be integrated over the relevant matrix ensemble. So solving a matrix model 
really amounts to computing integrals over matrix ensembles. 

The lectures will be divided into two steps: first we show how to interpret such matrix 
integrals in terms of discrete two-dimensional quantum gravity, namely in terms of graphs 
with prescribed topology and valences, carrying also configurations of statistical "matter" 
models; in a second step, we show how to compute these integrals explicitly. The main 
difficulty here is that the immense power of matrix integrals allows to get right and simple 
answers, but gives no really good reason for such simplicity, except for technical miracles 
that are sometimes called "integrability" . To compensate for this lack of understanding, 
we will always try to develop parallelly to the matrix model techniques and calculations 
some purely combinatorial reading of the various results. 

The simplest combinatorial objects in many respects are trees, and we will see, at least 
in the planar case, how graphs representing discrete surfaces of genus zero are reducible 
to decorated trees. This eventually explains the simplicity of the corresponding matrix 
model results. By pushing these ideas a little further, we will be able to investigate refined 
properties of discrete surfaces (graphs), involving their intrinsic geometry. For instance we 
will compute correlation functions for surfaces with marked points at a prescribed geodesic 
distance from one-another. 

Having collected many exact solutions for models of discrete geometry, it is natural to 
go to the continuum limit, which displays a rich singularity structure: indeed singularities 
may arise from the graphs themselves, say when parameters coupled to valences reach some 



critical values, and the contribution from large graphs start dominating the statistical sum. 
They may also arise from criticity of the matter statistical models defined on these already 
critical graphs, in which case collective behaviors start dominating configurations. The 
matrix models allow for taking both limits simultaneously (the so-called double-scaling 
limit) while keeping track of all genera. The continuum model is expected to be described 
by conformally invariant matter field theories [fE] coupled to 2D quantum gravity, i.e. 
defined on random surfaces || . Similarly we will write continuum correlation functions of 
the geodesic distance on the corresponding random surfaces. 

1.2. A brief history 

Planar graphs first arose in combinatorics, in the groundbreaking works of Tutte 0] in 
the 60's, who was able to compute generating functions for many classes of such objects, 
usually called maps by combinatorists. Higher genus was not considered then, and came up 
only later in physics works. The intrusion of matrix models in this subject occurred with 
the fundamental observation, due to t'Hooft || in the 70's, that planar graphs appearing 
in QCD with a large number of colors could be viewed as Feynman diagrams for matrix 
models, and that moreover the size of the matrices could serve as an expansion parameter 
to keep track of the topology of these diagrams. This caused the interest for matrix model 
to immediately rise, and led to the basic work of Brezin, Itzykson, Parisi and Zuber ||, 
who used various techniques to compute these matrix integrals, and among other things 
made the contact with Tutte's enumeration results. The matrix model techniques were 
then perfected by a number of people, whose list would be too long and probably not 
exhaustive. Then came the invention of continuum and discrete quantum gravity || , as the 
coupling of matter theories to fluctuations of the underlying space, both in field-theoretical 
and matrix languages. This second life of matrix models came to a climax in 1990 with 
three quasi-simultaneous papers making drastic progress in two-dimensional quantum 
gravity, as a toy model for low-dimensional non-critical strings, via the double-scaling limit 
of matrix models. This started a new matrix crazyness, and certainly helped develop matrix 
model theory a great deal (see || for a review and references). Remarkably, new areas of 
mathematics got infected by the matrix virus, thanks to Witten and Kontsevich @, who 
formulated a mathematically rigorous approach to the moduli space of punctured Riemann 
surfaces using matrix models, and set the ground for a little revolution in enumerative 
geometry. 



On the combinatorics front, it was only recently understood how to continue Tutte's 



work for higher genus graphs or more complicated planar cases [1C], but a good relation to 



matrix model results is still to be found. For planar graphs however, the simplicity of the 



matrix model solutions has finally been explained combinatorially by Schaeffer [ll|, who 



found various bijections between planar graphs and trees, allowing for a simple enumera- 



tion, and a precise contact with the matrix model solutions [fl2 ] . A remarkable by-product 



of this approach is that one may keep track on the trees of some features of the planar 



graphs, such as geodesic distances between vertices or faces []T3j [Tj], a task beyond the 
reach of matrix models so far. 



2. Matrix models for 2D quantum gravity 

2.1. Discrete 2D quantum gravity 

The purpose of quantum gravity is to incorporate in a field-theoretical setting the in- 
teractions between matter fields and the fluctuations of the underlying space. In Euclidian 
2D quantum gravity, the latter are represented by dynamical surfaces £ endowed with a 
Riemannian metric g and scalar curvature R, and for which the Einstein action of General 
Relativity reads 

S E = A [ yfgtft + M [ ^gRd 2 £ 

is is (2.1) 

= AA(E)+Af X (E) 

made of a cosmological term, in which the coslological constant A is coupled to the area of 
the surface A(E) and of the Newton term, in which the Newton constant M is coupled to 
the Euler characteristic x(S) of the surface. The dynamical surfaces are then discretized 
in the form of graphs with prescribed topology. 

We will now explain how matrix integrals can be used to generate such graphs, while 
precisely keeping track of their area and their Euler characteristic. For pedagogical pur- 
poses, we start with some simple remark on ordinary Gaussian integration, before going 
into the diagrammatics of Gaussian matrix integrals. 



2.2. Gaussian integral's diagrammatics 

Consider the following Gaussian average 
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e~^x 2n dx = (2n-l)! 



(2n)! 



(2-2) 



/2itJ-oc 2-n! 

Among the many ways to compute this integral, let us pick the so-called source integral 
method, namely define the source integral 
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(2.3) 



Then the average ( |2.2| ) is obtained by taking 2n derivatives of S(s) = e~ w.r.t. s and 
by setting s = in the end. It is then immediate to see that these derivatives must be 
taken by pairs, in which one derivative acts on the exponential and the other one on the 
prefactor s. Parallelly, we note that (2n — 1)!! = [2n — l){2n — 3). ..3.1 is the total number 
of distinct combinations of 2n objects into n pairs. We may therefore formulate pictorially 
the computation of (|2.2|) as follows. 



2n 1 





Fig.l: A star-diagram with one vertex and 2n out-coming half-edges stands 
for the integrand x 2n . In the second diagram, we have represented one non- 
zero contribution to (x 2n ) obtained by taking derivatives of S(s) by pairs 
represented as the corresponding pairings of half-edges into edges. 

We first draw a star-graph (see Fig.l), with one central vertex and 2n outcoming half- 
edges labelled 1 to 2n clockwise, one for each x in the integrand (this amounts to labelling 
the x's in x 2n from 1 to 2n). Now the pairs of derivatives taken on the source integral 
are in one-to-one correspondence with pairs of half-edges in the pictorial representation. 
Moreover, to get a non-zero contribution to (x 2n ), we must saturate the set of 2n legs by 
taking n pairs of them. Let us represent each such saturation by drawing the corresponding 
edges as in Fig.l. We get exactly (2n — 1)!! distinct labeled closed star-graphs with one 
vertex. This is summarized in the one-dimensional version of Wick's theorem: 



pairings 



(2.4) 



where the sum extends over all pairings saturating the 2n half-edges, and the weight is 
simply the product over all the edges thus formed of the corresponding averages (x 2 ) = 
(<i 2 /(is 2 )E(s)| s= o = 1. Each saturation forms a Feynman diagram of the Gaussian average. 
The edge pairings are propagators (with value 1 here) . This may appear like a complicated 
way of writing a rather trivial result, but it suits our purposes for generalization to matrix 
models and graphs. 

2.3. Gaussian matrix integral and more diagrammatics 

Let us now repeat the calculations of the previous section with the following Gaussian 
Hermitian matrix average of an arbitrary function / 

(f(M)) = ^y JdMe- NTr ^f(M) (2.5) 

where the integral extends over Hermitian NxN matrices, with the standard Haar measure 
dM = Y\ i dMaY\ i< -dRe{Mij)dIm{Mij), and the normalization factor Zq(N) is fixed by 
requiring that (1) = 1 for / = 1. Typically, we may take for / a monomial of the 
form f(M) = Ylfi j) e iMij, I a finite set of pairs of indices. Note the presence of the 
normalization factor N (=the size of the matrices) in the exponential. Note that the case 
of the previous section is simply the particular case of integration over lxl Hermitian 
matrices (i.e. real numbers) here. 

Like before, for a given Hermitian NxN matrix S, let us introduce the source integral 

E{S ) = ( e Tr(SM) )=e IrLp (2 _ 6) 

easily obtained by completing the square M 2 - N(SM + MS) = (M - NS) 2 - N 2 S 2 
and performing the change of variable M' = M — NS. We can use ( |2.6| ) to compute any 
average of the form 



<"«"«-> = aaraaP" E(s) l« (2 ' 7) 



d d_ 

3ji dSik 

Note the interchange of the indices due to the trace Tr(MS) = Yl^ijSji- As before, 
derivatives w.r.t. elements of S must go by pairs, one of which acts on the exponential and 
the other one on the S element thus created. In particular, a fact also obvious from the 
parity of the Gaussian, (277) vanishes unless there are an even number of matrix elements 



of M in the average. In the simplest case of two matrix elements, we have 



{M " Ma) = af^ SyeZ# 
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= irftiuSjk (2.1 

s=o JV 



Hence the pairs of derivatives must be taken with respect to Sij and Sji for some pair i, j 
of indices to yield a non-zero result. This leads naturally to the Matrix Wick's theorem: 



( n m ^ = e n ( m *j m *i) 

(i,j)El pairings P (ij),(kl)eP 



(2.9) 



where the sum extends over all pairings saturating the (pairs of) indices of M by pairs. 

We see that in general, due to the restrictions ( |2.8| ) many terms in ( |2.9| ) will vanish. 
Let us now give a pictorial interpretation for the non- vanishing contributions to ( |2.9| ) . We 
represent a matrix element Mjj as a half-edge (with a marked end) made of a double-line, 
each of which is oriented in an opposite direction. We decide that the line pointing from 
the mark carries the index z, while the other one, pointing to the mark, carries the index 
j. This reads 



M. 



Kl 



-<- 



(2.10) 



The two-element result (|2.8| ) becomes simply the construction of an edge (with both ends 
marked) out of two half-edges M{j and M^i , but is non-zero only if the indices i and j are 
conserved along the oriented lines. This gives pictorially 



(MijMji) 



m.. : > : !f T AT 



(2.11) 



Similarly, an expression of the form Tr(M n ) will be represented as a star-diagram with 
one vertex connected to n double half-edges in such a way as to respect the identification 
of the various running indices, namely 



Tr(ikT) = J2 M ili2 M i2i3 ...M inil 



«l,!2,-. ! n 




;2.i2) 





(a) 



(b) 



Fig. 2: An example of planar (petal) diagram (a) and a non-planar one (b). 
Both diagrams have n = 2p = 12 half-edges, connected with p = 6 edges. 
The diagram (a) has p + 1 = 7 faces bordered by oriented loops, whereas (b) 
only has 3 of them. The Euler characteristic reads 2 — 2h = F — E + 1 (V = 1 
in both cases), and gives the genus h = for (a), and h = 2 for (b). 

As a first application of this diagrammatic interpretation of the Wick theorem ( |2.9| ) , 
let us compute the large N asymptotics of (Tr(M n )). To compute (Tr(M n )), we must first 
draw a star-diagram as in ( |2.12| ), then apply ( |2.9|) to express the result as a sum over the 
saturations of the star with edges connecting its outcoming half-edges by pairs. To get a 
non-zero result, we must clearly have n even, say n = 2p. Again, there are (2p — 1)!! such 
pairings, and indeed we recover the case of previous section by taking N — 1. But if instead 
we take N to be large, we see that only a fraction of these (2p — 1)!! pairings will contribute 
at leading order. Indeed, assume first we restrict the set of pairings to planar ones (see 
Fig. 2 (a)), namely such that the saturated star diagrams have a petal structure in which 
the petals are either juxtaposed or included into one-another (with no edges-crossings). 
We may compute the genus of the petal diagrams by noting that they form a tessellation 
of the sphere (=plane plus point at infinity). This tessellation has V = 1 vertex (the star), 
E = p edges, and F faces, including the "external" face containing the point at infinity. 
The planarity of the diagram simply expresses that its genus h vanishes, namely 



2-2h = 2 = F-E + V = F + l-p =► F=p+l 



(2.13) 



Such diagrams receive a total contribution 1/N P from the propagators (weight 1/N per 
connecting edge), but we still have to sum over the remaining matrix indices ji,J2, ■■■■,jp+i 
running over the p + 1 oriented loops we have created, which form the boundaries of the 
F = p + 1 faces. This gives a weight iV per face of the diagram, hence a total contribution 
of iV p+1 . So all the petal diagrams contribute the same total factor N p+1 /N p = N to 



(Tr(M n )). Now any non-petal (i.e. non-planar, see Fig. 2 (b)) diagram must have at least 
two less oriented loops. Indeed, its Euler characteristic is negative or zero, hence it has 
F < E — V = p — 1 and it contributes at most for N F ~ P < 1/N. So, to leading order in N, 
only the genus zero (petal) diagrams contribute. We simply have to count them. This is 
a standard problem in combinatorics: one may for instance derive a recursion relation for 
the number c p of petal diagrams with 2p half-edges, by fixing the left end of an edge (say at 
position 1), and summing over the positions of its right end (at positions 2j, j = l,2,...,p), 
and noting that the petal thus formed may contain Cj—i distinct petal diagrams and be 
next to c p -j distinct ones. This gives the recursion relation 

p 
c p = ^2 c j-i c p-j c o = 1 (2.14) 

i=i 

solved by the Catalan numbers 

c p = . v V XI , (2.15) 

Finally, we get the one-matrix planar Gaussian average by taking the large iV limit: 

lim — (Tr(M n )) = ( C p lf n = 2 P ( 2 .16) 

n^oo N [0 otherwise v ' 

This exercise shows us what we have gained by considering N x N matrices rather than 
numbers: we have now a way of discriminating between the various genera of the graphs 
contributing to Gaussian averages. This fact will be fully exploited in the next example. 

2.4- Model building I: using one-matrix integrals 



Let us apply the matrix Wick theorem ( |2.9| ) to the following generating function 
f(M) = exp(A^J^ i>1 (7jTr(M l )/i), to be understood as a formal power series of the gi, 



i = l,2,3,4,... 



Z N ( 9l , 92 ,...) = { e N ^>-^^) 



*-*i 1J - i ni ny. LL (o -> 7 \ 

y- TT ( iV ^) ni y- n -e(t) n f(t) 

71l,n2,...>0 i>l ' a11 labelled fatgraphs V 

with 71,' i — valent vertices 




Fig. 3: A typical connected fatgraph T, corresponding to the average 
(Tr(M) 3 Tr(M 2 ) 2 Tr(M 3 )Tr(M 4 ) 2 Tr(M 6 )Tr(M 8 )). The graph was obtained 
by saturating the ten star-diagrams corresponding to the ten trace terms, 
namely with ri\ = 3 univalent vertices, n^ = 2 bi-valent ones, n^ = 1 tri-valent 
one, ri4 = 2 four-valent ones, rig = 1 six-valent one and n& — 1 eight-valent 
one, hence a total of y = 10 vertices. This graph corresponds to some par- 
ticular Wick pairing for which we have drawn the E — 16 connecting edges, 
giving rise to F = 2 oriented loops bordering the faces of V. 



by direct application 

In (|2.17| ), we have first represented pictorially the integrand ^(T^M 1 )) 71 * as a succession 

of rii i-valent star diagrams like that of Q2.12 ), i = 1,2, .... Then we have summed over 



all possible saturations of all the marked half-edges of all these stars, thus forming (non- 
necessarily connected) ribbon or fatgraphs F with some labelling of their half-edges (see 
Fig. 3 for an example of connected fatgraph). In (|2.17|) , we have denoted by E(F) the total 



number of edges of T, connecting half-edges by pairs, i.e. the number of propagators needed 
(yielding a factor 1/N each, from ( |2.8| )). The number F(T) is the total number of faces of 
r. The faces of T are indeed well-defined because T is a fatgraph, i.e. with edges made of 
doubly oriented parallel lines carrying the corresponding matrix indices i = 1,2, ...N: the 
oriented loops we have created by the pairing process are interpreted as face boundaries, 
in one-to-one correspondence with faces of F. But the traces of the various powers of M 
still have to be taken, which means all the indices running from 1 to N have to be summed 
over all these loops. This results in the factor iV per face of T in ( |2.17| ). Finally, the sum 
extends over all (possibly disconnected) fatgraphs T with labelled half-edges. Each such 
labelled graph corresponds to exactly one Wick pairing of fl2.9| ). Summing over all the 
possible labellings of a given un-labelled fatgraph F results in some partial cancellation 
of the symmetry prefactors ]^ l/(i ni rii\), which actually leaves us with the inverse of the 

10 



order of the symmetry group of the un-labelled fatgraph T, denoted by l/|Au£(r)|. This 
gives the final form 

N V(T)-E(T)+F(r) 

Z N ( 91 ,92,...)= £ n^' (r) (2-18) 

fatgrapha ' \ / I j>^ 

r — 

where 7i»(r) denotes the total number of i-valent vertices of F and V(T) = J^. 7*t(r) is the 
total number of vertices of I\ To restrict the sum in ( |2.18| ) to only connected graphs, we 
simply have to formally expand the logarithm of Z^, resulting in the final identity 

N 2-2h(T) 



connected 
fatgraphs T 



\Aut(F)\ 



where we have identified the Euler characteristic x(T) = F — E + V = 2 — 2h(F), where 
h(F) is the genus of F (number of handles). Eqn.( |2"TT9|) gives a clear geometrical meaning 
to the Gaussian average of our choice of f(M): it amounts to computing the generating 
function for fatgraphs of given genus and given vertex valencies. Such a fatgraph F is in 
turn dual to a tessellation F* of a Riemann surface of same genus, by means of rii z-valent 
polygonal tiles, i = 1, 2, .... 



The result (|2.19|) is therefore a statistical sum over discretized random surfaces (the 
tessellations), that can be interpreted in physical terms as the free energy of a model 
of discrete 2D quantum gravity. It simply identifies the Gaussian matrix integral with 
integrand f(M) as a discrete sum over configurations of tessellated surfaces of arbitrary 
genera, weighted by some exponential factor. More precisely, imagine only gs = g ^ while 
all other g^s vanish. Then ( |2.19|) becomes a sum over fatgraphs with cubic (or 3-valent) 
vertices, dual to triangulations T of Riemann surfaces of arbitrary genera. Assuming these 
triangles have all unit area, then ns(F) = A(T) is simply the total area of the triangulation 
T. Hence (|2.19| ) becomes 

^ n A (T) N 2-2h(T) 

connected triangulations T 

and the summand g N 2 ~ 2h = e~ E is nothing but the exponential of the discrete version 
of Einstein's action for General Relativity in 2 dimensions ( |2.1| ), in which we have identified 
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(a) 



(b) 



■4 1 




4 1 





Fig. 4: A 4-valent planar graph with hard dimers, represented by thickened 
edges. The corresponding graph obtained by shrinking the dimers (b) has 
both 4-valent and 6-valent vertices. The correspondence is three-to-one per 
dimer, as shown. 



the two invariants of E: its area A(E) and its Euler characteristic x(^) = 2 — 2/i(E). The 
contact with (|2.20| ) is made by setting g = e~ A and N = e 



-M 



If we now include all g^s in ( |2.19| ) we simply get a more elaborate discretized model, 
in which we can keep track of the valencies of vertices of T (or tiles of the dual T*). These 
in turn may be understood as discrete models of matter coupled to 2D quantum gravity. 
This is best seen in the case of the Hard-Dimer model on random 4-valent graphs [i5f. The 
configurations of the model are made of arbitrary 4-valent fatgraphs of arbitrary genus (the 
underlying discrete fluctuating space) and of choices of edges occupied by dimers, with the 
hard-core condition that no two adjacent edges may be simultaneously occupied (see Fig. 4 
for an illustration in the case of a planar graph). These matter configurations are given an 
occupation energy weight z per dimer, while the space part receives the standard weight 
g per 4-valent vertex, and the overall weight N 2 ~ 2h for each graph of genus h. We then 
note that any occupied dimer may be shrunk to naught, thus creating a 6-valent vertex by 
the fusion of its two 4-valent adjacent vertices. Comparing the configurations of the Hard- 
Dimer model on 4-valent graphs and those of graphs with only 4- and 6-valent vertices, 
we see that there is a one-to-three correspondence between those, as there are exactly 
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three ways of decomposing a 6-valent vertex into two adjacent 4-valent ones connected by 
a dimer (see the bottom line of Fig. 4). The Hard-Dimer model is therefore generated by 
an integral of the form ( |2.18| ), with only g± and g§ non-zero, and more precisely g& = g 
and #6 = 3g 2 z (=three decompositions into two 4-valent vertices and one dimer). This is 
the simplest instance of matter coupled to 2D quantum gravity we could think of, and it 
indeed corresponds to graphs with specific valence weights. 

Going back to the purely mathematical interpretation of (2.19), we start to feel how 



simple matrix integrals can be used as tools for generating all sorts of graphs whose duals 
tessellate surfaces of arbitrary given topology. The size iV of the matrix relates to the 
genus, whereas the details of the integrand relate to the structure of vertices. An important 
remark is also that the large iV limit of ( |2.19| ) extracts the genus zero contribution, namely 
that of planar graphs. So as a by-product, it will be possible to extract results on planar 
graphs from asymptotics of matrix integrals for large size N. 

2.5. Model building II: using multi-matrix integrals 

The results of previous section can be easily generalized to multiple Gaussian integrals 
over several Hermitian matrices. More precisely, let Mi, M2, ... M p denote p Hermitian 
matrices of same size iV x N, and Q a ,b, a, 6 = 1, 2, ...,p the elements of a positive definite 
form Q. We consider the multiple Gaussian integrals of the form 

fdMi dM e~^~^> b =i TT{MaQabMb) f(Mi M) 
(/(Mr, ..., M p )) = J dM - dM ^ TrU / Q mY" (2 ' 21) 

The one-Hermitian matrix case of the previous section corresponds simply to p = 1 and 
Qi t i = 1. The averages ( |2.21|) are computed by extending the source integral method of 



previous section: for some Hermitian source matrices Si, ..., S p of size N x N, we define 
and compute the multi-source integral 

E(Si, ..., S p ) = (e^Li Tr ( s » M »)) = e ^ E", b=1 r rr(s a (Q- 1 )a, b s b ) ^^ 



and apply multiple derivatives w.r.t. to S a 's to compute any expression of the form ( |2.21|) 



before taking S a — > 0. As before, derivatives w.r.t. elements of the S's must go by pairs to 
yield a non-zero result. For instance, in the case of two matrix elements of M a 's we find 
the propagators 

((M a )y(M 6 )fe,) = ^MQ _1 )a,6 (2.23) 
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In general we will apply the multi-matrix Wick theorem 

( J] (M a ) tJ )= J2 II ((Mah(M b ) kl ) (2.24) 

(a i l')f J pairings pairs 

v ' P {aij),(bkl)eP 

expressing the multi-matrix Gaussian average of any product of matrix elements of the M's 
as a sum over all pairings saturating the matrix half-edges, weighted by the corresponding 
value of the propagator ( |2.23| ). Note that half-edges must still be connected according 
to the rule ( |2.11| ), but that in addition, depending on the form of Q, some matrices may 
not be allowed to connect to one another (e.g. if (Q _1 ) a & = for some a and 6, then 
(M a Mf,) = 0, and in such a case, there cannot be any edge connecting a matrix with index 
a to one with index b). 

This gives us much freedom in cooking up multi-matrix models to evaluate generating 
functions of graphs with specific decorations such as colorings, spin models, etc... This is 
expected to describe the coupling of matter systems (e.g. a spin model usually defined on 
a regular lattice) to 2D quantum gravity (by letting the lattice fluctuate into tessellations 



of arbitrary genera). Famous examples are the O(n) model [pT6 | , the q-states Potts model 



|17|| , both including the Ising model as particular cases. Other models of interest require 



to use different types of matrices, to best represent their degrees of freedom. This is the 
case for the 6 vertex model expressed in terms of complex matrices, and for the so-called 
IRF (interaction round a face) models, expressed in terms of complex rectangular arrays 

ra nn 



3. The one-matrix model I: large N limit and the enumeration of planar graphs 

In this section, we will mainly cover the one-matrix integrals defined in Sect. 2. 4. Multi- 
matrix techniques are very similar, and we will present them in a concluding section. More 
precisely, we will study the one-matrix integral 

. . JdMe-^VM . . 

Z ^ V ) = J dMe -NTrV (M) (3- 1 ) 

with an arbitrary polynomial potential, say 

n*) = y-Ef^> and %(*) = y (3.2) 

This contains as a limiting case the partition function ( |2.17| ) of Sect. 2. 4. Note also that we 
are not worrying at this point about convergence issues for these integrals, as they must 
be understood as formal tools allowing for computing well-defined coefficients in formal 
series expansions in the g's. 
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3.1. Eigenvalue reduction 

The step zero in computing the integral ( |3.1|) is the reduction to a iV-dimensional 
integral, namely over the real eigenvalues m\, ..., m^v of the Hermitian matrix M. This is 
done by performing the change of variables M — > (m, U), where m =diag(mi, ..., mjv), and 
U is a unitary diagonalization matrix such that M = UmU* , hence U G U(N)/U(1) N as 
U may be multiplied by an arbitrary matrix of phases. The Jacobian of the transformation 
is readily found to be the squared Vandermonde determinant 

J = A(m) 2 = ] J {mi - rrij) 2 (3.3) 

l<i<j<N 

A simple derivation consists in expressing the differential dM in terms of dU and dm in 
the vicinity of U = I, namely dM = dllm + dm + mdU^ , but noting that UU* = /, we 
get dU' — —dU, and finally dM = dm + [dU, m], or dMij = drriiSij + (mi — mj)dUij, from 
which we directly read the Jacobian ( |3.3| ). Performing the change of variables in both the 
numerator and denominator of (13.11) we obtain 



, N L.JV dmi...dm N A(m) 2 e 2-n=i ^ mi > 
Zn(V) = l* ^-LJ -—s- (3.4) 



J uN dmi...dmNA{fn) 2 e 2-~n=\ 
3.2. Large size: the saddle-point technique 

Starting from the iV-dimensional integral (|3.4|) , we rewrite 



, . f dm 1 ...dm N e~ N S(m 1 ,...,m N ) 

Z n(V) = J rj , M2q( r 3.5 

j dm 1 ...dm N e~ N s o(m 1: ...,m N ) 



where we have introduced the actions 

1 N 1 

S{m 1 ,...,m N ) = —22v(m i )-j^ Yl Log|mi-ra./| 

1 N 1 

Sb(mi,...,mjv) = jyz2 v o(. m i) ~ jp 5Z Log\mi-mj\ 



(3.6) 



For large N the numerator and denominator of Q3.5| ) are dominated by the semi-classical 
(or saddle-point) minima of S and So respectively. For S, the saddle-point equations read 

9S =* V\m 3 ) = 1 £ - J — (3-7) 

\ < * rn ■ — rn ■ 



dm-: N ^— ' m-j — m 

J l<i<N J 
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for j = 1, 2, ..., N. Introducing the discrete resolvent 

1 N 1 



A" ^^ z — mi 

evaluated at the solution mi,..,mjv to fl3~77[) , multiplying (|3.7|) by l/(N(z — rrij)) and 
summing over j, we easily get the equation 

A *-^ z — rrij 
j=i J 

= -L y I (_L L_) 

N 2 ^-— ' rrij — rrij \ z — rrij z — m* J / n n \ 

i<i^j<N J l v J l/ (3.9) 

= i y I 

A^ 2 ^— ' (z — mA(z — rrij) 



= uj n (z) 2 + —uj' n (z) 

Assuming uj^ tends to a differentiable function u(z) when JV ^ oo we may neglect the 
last derivative term, and we are left with the quadratic equation 

u(z) 2 -V'{z)uj{z)+P{z) = 

PW-ito, If "'(«>- "'<"*> (3 ' 10) 

N^oo N ^ z -rrij 

3=1 

where P(z) is a polynomial of degree d — 2,d the degree of V. The existence of the limiting 
resolvent u(z) boils down to that of the limiting density of distribution of eigenvalues 



1 N 
p(z) = lim — 22 6 ( z ~ m j) ( 3 - n ) 



normalized by the condition 



p(z)dz = 1 (3.12) 

as there are exactly A" eigenvalues on the real axis. This density is related to the resolvent 
through 

u(z) = f ^-dx = J2 — I x m - 1 p(x)dx (3.13) 
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where the expansion holds in the large z limit, and the integral extends over the support 
of p, included in the real line. Conversely, the density is obtained from the resolvent by 
use of the discontinuity equation across its real support 

p(z) = ■—— lim u>(z + ie) — u(z — ie) z G supp(p) (3-14) 

2l7T e^O 

Solving the quadratic equation ( |3.10|) as 

wW = n^WTO-4Pw (3 15) 

we must impose the large z behavior inherited from ( p . 1 2| ) ( |37l3|) , namely that u)(z) ~ 1/z 
for large z. For d > 2, the polynomial in the square root has degree 2{d — 1): expanding 
the square root for large z up to order 1/z, all the terms cancel up to order with V'(z), 
and moreover the coefficient in front of 1/z must be 1 (this fixes the leading coefficient of 
P). The other coefficients of P are fixed by the higher moments of the measure p(x)dx. 
For instance, when k = 2 and V — Vq, we get P — 1 and 

Wo (2) = ^-^ 2 -4) (3.16) 

It then follows from ( |3.14j) that the density has the compact support [—2, 2] and has the 
celebrated "Wigner's semi-circle law" form 

Po(z) = ^-a/4^ (3.17) 

The resolvent ujq is the generating function for the moments of the measure whose density 
is po ( y i a the expansion ( ^.13| )), from which we immediately identify 

*•*(*)*={* n:z (3.i 8 ) 



with c p as in ( |2.15| ). Indeed, due to the quadratic recursion relation ( |2.14| ), the gen 



erating function C(x) = J2 p >o xPc p satisfies xC(x) 2 = C(x) — 1, and therefore we have 
uiq{z) = C(l/z 2 )/z. The coefficients ( 3.18 ) are nothing but the planar limit of the Gaussian 
Hermitian matrix averages (with potential Vq(x) = x 2 /2), namely liniAr^ 00 (-^TrM n )vo = 
J R x n po(x)dx, hence our analytical result (|3.18|) is an alternative for that already obtained 
combinatorially in ( |2.16|) . 
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In the general case, the density reads 



1 



p(z) = ^V4P(z)-(V'(z)r (3.19) 

and may have a disconnected support, made of a union of intervals (the so-called multicut 
solutions). It is however interesting to restrict oneself to the case when the support of 
p is made of a single real interval [a, b], as this will always be the preferred saddle-point 
solution for generating the correct formal series expansions of the all-genus free energy. 
For supports made of more than one interval, resonances may occur as eigenvalues tunnel 
from one interval to another, and oscillations develop in the N dependence, which cause 
the large N expansion to break down, unless some strong conditions are imposed on say 
complex contour integrals for the eigenvalues. The one-cut hypothesis will be justified 
a posteriori in Sect. 4 below, when we revisit the problem from a purely combinatorial 
perspective. 

In the one-cut case, the polynomial V'(z) 2 — 4P(z) has single roots at say z = a 
and z = b and all other roots have even multiplicities. In other words, we may write the 
limiting resolvent as 

u(z) = \{V'{z) - Q(z)y/{z-a){z-b)) (3.20) 

where Q(z) is a polynomial of degree k — 2, entirely fixed in terms of V by the asymptotics 
oj(z) ~ 1/z for large \z\. More precisely, let us introduce H(z) = V' (z) / y/Jz — a)(z — 6), 
considered as a series expansion for large z, then Q(z) is nothing but the part of this series 
that is polynomial in z, denoted as H + (z). Writing moreover H(z) = H + (z) + H_(z), we 
finally get 

uj(z) = ^H_ (z) y/( z -a)(z-b) (3.21) 

Writing H-(z) = E l >i H- iZ ~\ we get that u(z) ~ 1/z iff H- X = and #_ 2 = 2. These 
coefficients are expressed as residue integrals at infinity, namely 

H_ m ( z ) = I ^ z m-i^^ M (3.22) 

The square root term is uniformized by the change of variables z = w + S + R/w, with 
S = z±!> and R = (^) 2 , and 

H- m {z) = I -^-{w + S + R/wr^V'iw + S + R/w) (3.23) 

J Z171W 
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so that finally if_i = Vq and H_ 2 = V!_ 1 J r SVq + RV{, where the shorthand notation V^ 
stands for the coefficient of w m in the large w expansion of V (w + S + r/w). Performing 
the change of variables w — > R/w allows to relate V'_ m = R m V^ n . Finally, the asymptotic 
condition w(z) = 1/z + 0(l/z 2 ) at large z boils down to 



K*-i)/2] ,, 1){ 



V' = = S-Y 9i y S'- 2 ^ 1 R \ — ¥ h rr 



i>\ j=0 

[i/2] 



(3.24) 



v 1 =i = R-y gi y &-**& .„ ( '~ 1)! — rr 



These equations simplify drastically in the case of even potentials, where gt = for all odd 
i. The parity of V indeed induces that of p, and we have S = (a + b)/2 = as the support 
of the density is symmetric w.r.t. the origin. This leaves us with only one equation 

l = R-y^g 2l R^ 2l : l \ (3.25) 

for R = a 2 /A. In the particular case of the gaussian potential V = Vq, this reduces to R = 1 
and S = 0, in agreement with b = —a = 2 ( |3.16| ). Expanding the solutions of ( |3.24[) as 
formal power series of the g^s, the conditions R = l + 0({gi}) and S = 0({gi}) determine 
them uniquely. These in turn determine a and b and therefore p and u completely. 

The planar free energy / = F - F = limjv^oo ^jLog (Z N (V)/Z N (V Q )) is finally 
obtained by substituting the limiting densities p, po i n the saddle point actions S and So, 
with the result F — Fq = Sq — S. It is however much simpler to evaluate some derivatives 
of the free energy, by directly relating them to the planar resolvent uj(z), the subject of 
next section. 

3. 3. Enumeration of planar graphs with external legs 

Let us first consider the generating function Ti for planar graphs with weights gi per 
i-valent vertex, and with one external (univalent) leg, represented in the external face on 
the plane (see Fig. 5 (a)): 

rx = df/dg x = lim l(Tr(M)) y = u. 2 = \(H_ Z - SH. 2 - 2RH- 1 ) (3.26) 



where use has been made of ( |3.13| ), and as before uj- m denotes the coefficient of z m in the 



large z expansion of uj (z). From the large z asymptotics of uj(z), we know that H-\ = and 
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(a) 



(b) 





(c) 



(d) 



Fig. 5: Samples of planar graphs with external legs (univalent vertices marked 
with a cross) and arbitrary valences, with respectively (a) one leg in the 
external face (b) one leg (anywhere) (c) two legs in the same (external) face 
(d) two- legs (one in the external face, the other anywhere). 

H-2 = 2, and we must now evaluate H_ 3 = VL 2 + 2SVL 1 + (S 2 + 2R)V^ + 2RSV{ + R 2 V 2 ' 
2{V'_ 2 + 2SVL 1 ) + (S 2 + 2R)V^ = 2(VL 2 + 2S), leaving us with 



r x = v'+s 



(3.27) 



Analogously, we may compute the connected two-leg-in-the-same-face graph generating 
function 1?2 = 1V-3 — T 2 (see Fig. 5 (c)), in which we subtract the contributions from 
disconnected pairs of one-leg graphs. We get uj-3 = df/dg 2 = R + S 2 + V!_ 3 + 2SV!_ 2 and 
finally 



r 2 = r + v!_ 3 
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(VL, 



(3.28) 



Another quantity of interest is the connected two-leg graph generating function T^i = 
d 2 f/dgf = duj-2/dgi (see Fig. 5 (d)). This turns into 

r,, = ^ « + ^» «(l + V?,) + f«V« (3.29) 

dgi dgi dg x dg x dg x 

Let us first replace the term 1 in factor of dS/dgi by 1 = V!_ x , the second equation of 
( p.24j) . Note that the residue of a total differential always vanishes, hence in particular 

(d/dw(wV'(w + S + R/w)))_ x =0 = 1/^ + V'l 2 - RV^ (3.30) 

This allows to rewrite 

r " = a|^" + ^ < 3 ' 31 ' 

Finally, differentiating the equation Vq = w.r.t. g x yields = dS/dgiVQ +dR/dgiV" — 1, 
where the last term comes from the explicit derivation w.r.t. g± of V'(x) = x — g\ — 
g 2 x — gsx 2 /2 — ... Multiplying this by R, and noting as before that RV{' = V" x , we get 
RdS/dg^' + dRjdgxV'lx = R and finally 



1,1 



R (3.32) 



This result holds for even potentials as well, upon setting all #2i+i = in the end. Eq. fl3.32p 
gives a straightforward combinatorial interpretation of R as the generating function for 
planar graphs with two external (univalent) legs, not necessarily in the same face. 

To conclude the section, let us now give a combinatorial interpretation for S. Let us 
show that S is the generating function for one-leg planar graphs. By this we mean that 
the leg need not be adjacent to the external face, as was the case for T x (see Fig. 5 (b)). 
Comparing with the definition of Ti, we must show that S is the generating function for 
one-leg planar graphs (with the leg in the external face), and with a marked face (chosen 
to be the new external face). This amounts to the identity 

S = zd z T 1 \ z=1 (3.33) 

where we have included a weight z per face of the graph, to be set to 1 in the end. Due 
to Euler's relation F = 2 + E — V, where E is the total number of edges, and V that of 
vertices of the one- leg graphs at hand, and noting that 2E = l + $^£Vi while V = 1 + ^2 V{, 
where V{ is the number of internal z-valent vertices, so that 2E — V = J^(i — 1)V^, we 
see that zd z Y\ = (2 + tdt)Ti, if we attach a weight 1/t per edge and t l_1 per i-valent 
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vertex (with a net resulting weight t 2E ~ v ~ E = t E ~ v ). Modifying the propagator and 
vertex weights of the matrix model accordingly, this simply amounts to replacing V'{x) by 
V'(tx) = tx — ^2 gj i £ l ~ 1 x l ~ 1 in all the above formulas, and setting t = 1 after differentiation. 
This yields 

aq ar> 

(2 + d t )r 1 \ t=1 = 2S + 2VL 2 + ^\t=ia + V!! 2 ) + ^\t=iV!! 1 + V!! 3 + SV!! 2 +RV!! 1 (3.34) 



We now use the above trick (]3.30|) that the residue of a derivative vanishes, but this time 
with 

(d/dw(w 2 V'(w + S + R/w)))^ = = 2V!_ 2 + V" 3 - RV!i x (3.35) 

and we use this to eliminate V" 3 from ( |3.34|) , as well as ( |3.30| ) to rewrite the factor of 
dS/dt as 1 + V!! 2 =VL X + V'l 2 = RV^', with the result 

aq ar> 

(2 + t )r!| t=1 = 2S + SV!! 2 + 2RV , : 1 + RV£'-^\t=i + V -i^\t=i (3-36) 

Let us now differentiate w.r.t. t the equation = Vq, and then set t = 1 and multiply it 
by R. This gives 



(3.37) 



aq flp 

= R(VH t + SV{' + RV X " + v^\ t=1 + v;'^\ t =i) 

aq ar> 

= RSV>> + 2RVH X + RV^\t=x + V-i^fU=i 
and allows to rewrite ( ^.36| ) as 

(2 + d t )T 1 \ t=1 = 2S + SVH 2 - RSVf = S + S(V!_ X + V'l 2 - RV{J) = S (3.38) 

by replacing 1 — ► V!_ x and using again the equation ( |3.30|) . This completes the identification 
of S as the generating function for one-leg planar graphs, with the leg not necessarily in 
the external face. 

That the generating functions for both one- and two-leg planar graphs should satisfy 
a system of two algebraic equations ( |3.24| ), looks like magic at first sight. It is the purpose 
of Sect. 4 below to unearth the combinatorial grounds for this apparent miracle. 
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3. 4- The case of 4-valent planar graphs 

Before going into this, let us conclude with the case of the quartic potential say 

O A 

V(z) 



4j — g^, for which we have S = and eq. (|3.25| ) reduces to 



1 = R - 3gR" 



R 



a 

T 



l 

6g' 



(1 - V 1 " 12 9) 



(3.39) 



as R is the unique solution with the power series expansion R = 1 + 0(g). The corre- 
sponding resolvent and density of eigenvalues read respectively 



u>(z) = ~(z 



gz 6 - (i 

,2 



#y -gz 2 )y/z 2 - 



o/ 



1 f z 

p(z) = — (1 - g— - gz 2 )Va 2 -z 2 

Ztv Z 



(3.40) 



The two- leg-in-the-same- face graph generating function 1?2 of eq.( |3~28|) reads here 

v 2 = R - gR ^ R ^; R ) 



(3.41) 



where we have used eq.( |3.39| ) to eliminate g. But any planar 4-valent graph with two 
external legs in the same face is obtained by cutting an arbitrary edge in any closed planar 
4-valent graph. As the two legs are distinguished, and as there are exactly twice as many 
edges than vertices in a closed 4-valent graph, we have 1?2 = l+4gdf /dg. The contribution 
1 comes from the unique graph made of one loop, with one edge and no vertex, not counted 
in /. This gives the differential equation 

df (R-l)(3-R) 



4g- 



(3.42) 



dg 3 

and eliminating g = (R-1)/(3R 2 ) from ( ^39| ), we finally get §- = R 2 (?>-R)/4. Changing 



dg 

variables to R, this turns into -Jb = (2 — -R)(3 — R)/(12R), easily integrated into 



an 



hogR+^(R-l)(R-9) 



(3.43) 



where the constant of integration is fixed by requiring that / = when R = 1 (Gaussian 
case V = Vq). Substituting the expansion R = 1 + 3g + 18g 2 + ... into ( |3.43|) yields the 
expansion 

+ ... (3.44) 



= g 



2 2 2 8 

where we have represented the planar 4-valent graphs with up to 2 vertices, together with 
their inverse symmetry factors. 
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4. The trees behind the graphs 

Using the above interpretation of R as the generating function for planar graphs with 
two distinguished external legs not necessarily in the same face, let us now establish a 
general bijection between such graphs and suitably decorated trees, also called blossom- 
trees. 

4-1. \-valent planar graphs and blossom trees 





(a) 



(b) 



±+' 




4^ 




(c) 



(d) 



Fig. 6: Illustration of the bijection between two-leg planar 4-valent graphs and 
rooted blossom trees. Starting from a two- leg graph (a), we apply the iterative 
cutting procedure, which here requires turning twice around the graph. In 
(b), the indices indicate the order in which the edges are cut during the 1st 
turn (1, 2, 3) and 2nd turn (4, 5, 6). Each cut edge is replaced by a black/white 
leaf pair (c) , while the in-coming leg is replaced by a leaf and the out-coming 
one by a root, finally leading to a blossom tree (d). Conversely, the matching 
of black and white leaves of the blossom tree (d) rebuilds the edges of (a). 
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For reasons of simplicity, let us start with the case of 4-valent graphs. Given a two- leg 
such graph G (see Fig. 6 for an illustration), we represent it in the plane by picking the 
external face to be adjacent to the first (in-coming) leg. We now visit all edges bordering 
this external face in counterclockwise direction, and cut them iff the resulting graph remains 
connected. We then replace the two halves of the cut edges by respectively a black and a 
white leaf. This "first passage" has merged a number of faces of the initial graph with the 
external one. We now repeat the algorithm with the new external face, and so on until all 
faces are merged. The resulting graph is a 4-valent tree T (by construction, it has only 
one face and is connected). The tree is then rooted at its second (outcoming) leg, while 
its incoming one is replaced with a white leaf. Attaching a charge +1 (resp. —1) to white 
(resp. black) leaves, we obtain a tree with total charge +1. It is easy to convince oneself 
that the resulting 4-valent tree has exactly one black leaf at each vertex. 



///// 




m 



m=p+l 




m 



Fig. 7: The only two possibilities for the environment of an edge in a 4-valent 
blossom-tree, obtained by cutting a two-leg planar 4-valent graph. The edge 
separates the tree into a top and a bottom piece. The first leg of the graph is 
chosen to be in the bottom piece. The two cases correspond to whether the 
cutting process stops in the top (a) or bottom (b) piece. We have represented 
in both cases only the leaves unmatched within each piece. In each case, the 
position of the root (second leg) is fixed by the fact that any 4-valent tree 
must have an even number of leaves (including the root). We have indicated 
the corresponding charges q = or +1 of the top and bottom pieces. 
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This is best proved by showing that its descendent subtrees not reduced to a black leaf all 
have charge +1. To see why, consider any edge of the blossom tree, not directly attached 
to a black leaf. It separates the tree into two (top and bottom) pieces as depicted in Fig. 7. 
As a result of the above iterative cutting procedure, we may keep track of the m and p 
cut edges encompassing this edge, respectively lying on its right and left, and connecting 
the top and bottom pieces. Assuming the first leg was in the bottom part, and as the 
cutting process travels in counterclockwise direction, we may only have m = p + 1 or 
m = p according to whether the cutting process stopped in the top or bottom piece. But 
as the top and bottom pieces are trees with only 4-valent inner vertices, they must have an 
even number of leaves, including the root, and the cut edge. Eliminating those matched 
by black/white pairs within each piece, we are respectively left with: in case (a), 2p + 2 
leaves on top and 2p + 3 on the bottom, hence the root must be in the bottom; in case (b), 
2p + l leaves on top and 2p + 2 on the bottom, hence the root must be on top. Adding up 
the charges, we see that the descendent piece (not containing the root) always has charge 
q = +l. 

Let us now define rooted blossom-trees as rooted planar 4-valent trees with black and 
white leaves, a total charge +1, and exactly one black leaf at each vertex (or equivalently 
such that each subtree not reduced to a black leaf has charge +1). Then the rooted blossom- 
trees are in bijection with the two-leg 4-valent planar graphs. The inverse mapping goes 
as follows. Starting from a rooted blossom-tree T, we build a two-leg 4-valent planar 
graph by connecting in counterclockwise direction around the tree all pairs of black/white 
leaves immediately following one-another, and by repeating this until all black leaves are 
exhausted. This leaves us with one unmatched white leaf, which we replace by the first leg, 
while the root becomes the second leg. The order in which leaves are connected exactly 
matches the inverse of that of the above cutting procedure. This bijection now allows for 
a direct and simple counting of 2-leg 4-valent planar graphs, as we simply have to count 
rooted blossom-trees. Decomposing such trees according to the environment of the first 
vertex attached to their root, we get the following equation for their generating function 



R 



+ 



3gR 2 



(4.1) 



um 



© 



+ g 
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where the first term corresponds to no vertex (and a white leaf directly connected to the 
root), and the three others to a vertex with one black leaf and two descendent blossom- 
trees, each receiving a weight g for the decomposed vertex. Note that eq.(fO) trivially 



amounts to the first equation of (|3.39|) . We have therefore found a purely combinatorial 
re-derivation of the one-cut large N matrix model result for planar 4-valent graphs, which 
confirms its validity. 

4-2. Generalizations 

More generally, the above bijection may be adapted to two-leg planar graphs with 
arbitrary even vertex valences. Repeating the above cutting procedure on such a two-leg 
planar graph leaves us with a rooted tree with only even vertex valences, with black and 
white leaves, and a total charge +1, but with now exactly k — 1 black leaves attached 
to each of its 2/c-valent vertices. This is again a consequence of the equivalent property 
that any subtree not reduced to a black leaf has charge +1, a fact proved exactly in the 
same manner as before (actually, Fig. 7 is still valid for the case of arbitrary even valences). 
This suggests a straightforward generalization of rooted blossom-trees with arbitrary even 
vertex valences, with black and white leaves and such that any subtree not reduced to a 
black leaf has charge +1. The latter are again in bijection with the two-leg planar graphs 
with even valences, and are easily enumerated by considering the environment of the vertex 
attached to the root, with the result for the generating function R, including weights gik 
per 2/c-valent vertex: 

R=l + Y,9*k( 2k ~ 1 )R h (4-2) 

fc>i ^ ' 



where the first term corresponds as in G4.lT ) to the tree with no vertex, while the /c-th term 
in the sum corresponds to the ( C" ) ways of picking the fe — 1 black leaves among the 
2k — 1 descendents of the 2/c-valent vertex attached to the root, the remaining descendents 
being themselves trees of charge +1 generated by R. The equation fl4.2| ) is nothing but 
( p.25| ), written in a different fashion. 

Finally, the bijection may be adapted so as to also include arbitrary (both even or 
odd) valences, but then requires the introduction of one-leg graphs as well. Such graphs are 
represented in the plane with their unique leg not necessarily adjacent to the external face, 
hence are not generated by Ti = S + V!_ 2 , but, as we showed in the previous section, by S 
itself. The graphs are again cut according to the above procedure, to produce rooted trees. 
The system of equations ( |3.24| ) is nothing but that obeyed by the rooted blossom-trees of 
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two kinds corresponding to cutting one- and two-leg graphs, respectively generated by S 
and R, and denned as rooted trees with black and white leaves, and total charge and +1 
respectively, and whose descendent subtrees not reduced to a black leaf all have charge 
or +1. A simple way of recovering all combinatorial factors in the two lines of ( |3.24|) is to 
note that in a rooted blossom tree of charge (resp. +1), the i — 1 descendents subtrees of 
any i-valent vertex attached to the root may be either black leaves (charge —1), blossom 
trees of charge 0, or blossom trees of charge 1, the total charge being (resp. +1). These 
subtrees are generated respectively by the functions 1, S and R. Denoting by j in both 
cases the total number of descendent subtrees of charge +1, we must have j (resp. j — 1) 
black leaves to ensure the correct total charge, and the remaining i — 2j — 1 (resp. i — 2j) 
descendents have charge 0. The combinatorial factors of ( |3.24| ) account for the possible 
choices of these among the i — 1 descendents. 

This combinatorial interpretation sheds light on the algebraicity of the equations ob- 
tained in the large iV limit for the general one-matrix model: trees are indeed archetypical 
objects whose generating functions obey algebraic relations, and we have shown that the 
planar graphs generated by the large iV matrix model could be represented by (blossom) 
trees. This correspondence will be fully exploited in Sect. 6 to investigate the intrinsic 
geometry of planar graphs. 

5. The one-matrix model II: topological expansions and quantum gravity 

We now turn to higher genus contributions to the one-matrix model free energy. This 
is best done by use of the so-called orthogonal polynomial technique [|20 . 



5.1. Orthogonal polynomials 



The standard technique of computation of ( |3.4j ) uses orthogonal polynomials. The idea 



is to disentangle the Vandermonde determinant squared interaction between the eigenval- 

m — 1 

j=0 P™,3 : 



ues. The solution is based on the following simple lemma: if p m (x) = x m + J2T=o Pm,j% 



are monic polynomials of degree m, for m — 0, 1, ..., N — 1, then 

A(ra) = det(ra^~ )i<i,j<N = det(pj_i(mi))i<jj<jv (5.1) 

easily derived by performing suitable linear combinations of columns. Let us now introduce 
the unique set of monic polynomials p m , of degree m = 0,1, ...,N — 1, that are orthogonal 
w.r.t. the real one-dimensional measure dfi(x) = exp(—NV(x))dx, namely such that 

(p m ,Pn)= / p m (x)p n (x)dfj,(x) = h m 5 mjn (5.2) 
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These allow us to rewrite the numerator of ( |3.4| ), using (5J.), as 



N r N-l 

J2 <vt)H d^m i )p a{l) _ 1 {m i )p T(l) . 1 {m i )e- NV{m ^=N\\[h (5.3) 

a,T<=S N »=1 R 3=0 

We may apply the same recipee to compute the denominator, with the result N\ n ? =o hj > 
where the hS are the squared norms of the orthogonal polynomials w.r.t. the Gaussian 



measure d[io(x) = exp(— Nx 2 /2)dx. Hence the /i's determine Zn(V) entirely through 

AT-l 



z »<y> = n tw < 5 - 4 > 

i=0 n i 

To further compute the /i's, let us introduce the two following operators Q an d -P, 
acting on the polynomials p m : 

QPm(x) = Xp m (x) 

n r , d ,, (5 ' 5) 

PPm{x) = -j-Pm{X) 

with the obvious commutation relation 

[P,Q] = 1 (5.6) 

Using the self-adjointness of Q w.r.t. the scalar product (/, g) = J f(x)g(x)dfj,(x), it is easy 
to prove that 

QPm(x) = Xp m (x) = p m+ l(x) + S m p m {x) + r m Pm-l{x) (5.7) 

for some constants r m and s m , and that s m = if the potential V(x) is even. The same 
reasoning yields 

r m - h m , m = 1,2, ... (5.8) 

"m— 1 

and we also set ro = /to fo r convenience. 

Moreover, expressing both (Pp m ,Pm) and (Pp m ,p m _i) in two ways, using integration 
by parts, we easily get the master equations 

m (V'(Q)p m ,p m - 1 ) 



N (pm-l,Pm-l) (5.9) 

= (V'(Q)p mi p m ) 
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which amount to a recursive system for s m and r m . Note that the second line of ( |5.9|) is 
automatically satisfied if V is even: it vanishes as the integral over IR of an odd function. 
Assuming for simplicity that V is even, the first equation of ( |5.9[ ) gives a non-linear re- 
cursion relation for the r's, while the second is a tautology, due to the vanishing of all the 

s's: 

m (V'(Q)p m ,Pm-i) v^ (Q 2k ~ 1 Pm,Pm-i) 

fc>l 

2k-2 (5.10) 



N (Pm-l,Pm-l) ^ (Pm-bPm-l) 

2fc-2 



fc>l paths p|p(l)=m, p(2fc — l)=m— 1 j=l 

p(i+l)-p(i)=±l 

where the sum extends over the paths p on the non-negative integer line, with 2k — 1 steps 
±1, starting at p(l) = m and ending at p(2k — 1) = m — 1, and the weight reads w(p, q) = 1 
if g = p + 1, and w(p,q) = r p if q = p — 1. For up to 6-valent graphs this reads 

— = r„(l - #2) - #4Tn(rn+l + r n + r n-l) 

- 5r 6 (r n+ ir n+2 + r n+ ir n _i + r n _ir n _ 2 + r 2 n + r 2 n+1 + r\_ x + 2r n (r n+1 + r n -i) 

(5.11) 

In general, the degree d of V fixes the number d — 1 of terms in the recursion. So, we 

need to feed the d — 2 initial values of ro, ri, r%, ..., rd-3 into the recursion relation, and we 

obtain the exact value of Z^iV) by substituting hi = r§r\...Ti in both the numerator and 

the denominator of (|5.4f) . Note that for Vq(x) = x 2 /2 the recursion (|5.9|) reduces simply 

to 

m _ {Qp { l\pt ] -i) (o) , , 

N~ (JO) jo) x " r - ^^ 



and therefore ftm = ^0 m^/N 771 = y2ixrn\ / N mJrl ' 2 . The p™ are simply the (suitably 
normalized) Hermite polynomials. 

Finally, the full free energy of the model (|3J]) reads 



N — 1 

N %-;.„ ._ Nr, 




F N (V) =LogZ N (V) = NLogr J — + ^ (N - i)Log— - (5.13) 



2tt 

i=l 



in terms of the r's. 
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5.2. Large N limit revisited 

In view of the expression ( |5.13 ), it is straightforward to get large iV asymptotics for 
the free energy, by first noting that as ho ~ J j, the first term in ( |5.13| ) doesn't contribute 
to the leading order A^ 2 and then by approximating the sum by an integral of the form 

i Ar_1 ' f 1 ( ) 

f = „^ jf E (i - >°^ = I *d - *> WV ( 5 -») 

where we have assumed that the sequence r* tends to a function r* = r{i/N) of the variable 
z = i/N when N becomes large. This assumption, wrong in general, basically amounts 
to the one-cut hypothesis encountered in Sect. 3. 2. The limiting function r(z) in ( |5.14j ) is 
then determined by the equations (|5.9| ), that become polynomial in this limit. In the case 
V even for instance, we simply get 

z = r{z) - J2 ^1 J92kr(z) k (5.15) 

The function r(z) is the unique root of this polynomial equation that tends to z for small 
z (it can be expressed using the Lagrange inversion method for instance, as a formal power 
series of the g's), and the free energy follows from ( |5.14j ). To relate this expression to our 
former results, let us again attach an extra weight t per face of the graphs. As before, it 
amounts to replacing V'(x) — > V'(tx) = tx — ^ fc>1 g2kt 2k ~ 1 ( ^ )x 2k ~ 1 , and to rescale 
/ — > t 2 f. Setting p(z) = t 2 r(z), we arrive at 

tz = p{z) - J2 ( 2 \~ 1 )92kp(z) k = <p(p(z)) (5.16) 

fc>i ^ ' 

and / = L dz(l — z)Log^^-, p(z) being determined by cp(p(z)) = tz. Let us perform in 
this integral the change of variables z — > p, with dz = <p'(p)/tdp, and integration bounds 
p(0) = and p(l) = r, solution of t — <^(r). We obtain: 

t 2 f= f 'dp? '(p)(t - <p(p)) Log-^ (5.17) 

Jo <PW 

We now take derivatives w.r.t. t: as the dependence on t is either via r or explicit in the 

integrand, there are only two terms involved. But the integrand vanishes at the upper 

bound, as t — ip{r) = 0, hence only the explicit derivative contributes, and we have 



d t (t 2 f) = [ dpip'(p) Log 
Jo 



VV> (5.18) 

d 2 (t 2 f) = d t vy'(r) Log -^- = Log - 
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Note that r may be interpreted in the light of Sect. 4 as the generating function for 
rooted blossom trees with a weight t per white leaf (easily read off the relation r = 
* + S/c>i #2fc( k ) rk )- F ma hy, setting t — 1, we may rewrite 

d 2 t (t 2 f)\ t=1 = Logi? = -Log(l - Y,92k( 2k " 1 )^" 1 ) (5.19) 

fc>i ^ ' 

as r reduces to R at t = 1. This expresses the generating function for planar graphs 
with even valences and with two distinct marked faces (as each derivative amounts to a 
marking) as the logarithm of the generating function for blossom trees. This formula will 
become combinatorially clear in Sect. 7. 2 below. 

5.3. Singularity structure and critical behavior 

In the even potential case, according to ( |5.18| ) , the singularities of r govern those of 
the free energy, r attains a first critical singularity at some t = t c where r = r c with 
(p(r c ) = t c and f'(r c ) = 0. We may then Taylor-expand 

t c -t = <p(r c ) - <p(r) =~\{v- r c ) We) + 0((v - r c ) 3 ) (5.20) 

As t is an activity per face of the graphs, we may consider the number of faces as a 
measure of the area of the associated discrete surface, therefore the singularity r s i ng ~ 
(t c — t) 1 ' 2 is immediately translated via ( |5.18|) into a singularity of the planar free energy 
fsing ~ (t c — ^) 2_7 , with a string susceptibility exponent 7 = —1/2. Alternatively, upon 
Laplace-transforming the result, this exponent also governs the large area behavior of 
Ja ~ const, t"" 4 /^ 3-7 , the planar free energy for fixed area (A=number of faces here). 
This is the generic singularity expected from a model describing space without matter, 
such as that of the pure 4-valent graphs studied above. 

We may reach more interesting multicritical points with different universality classes 
and exponents by fine-tuning the parameters g^k so as to ensure that a higher order 
singularity is attained at some t = t c such that r = r c , while <p'(r c ) = <p"(r c ) = ... = 
(p( m >(r c ) = 0, while (p( m+1 '(r c ) 7^ 0. Taylor-expanding now yields 

t c -t = v?(r c ) - vp(r) = - V {m + [y (r ~ ^) m+1 + 0((r - r c ) m + 2 ) (5.21) 

This translates into a singularity of the free energy with string susceptibility exponent 
7 = ~ m +i - This is characteristic of non- unitary matter conformal field theory with central 
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charge c(2, 2m + 1) coupled to 2D quantum gravity [[I ] 0. The first example of this is the 
Hard Dimer model introduced in Sect. 2. 4 above, for which 

^HD(r) = r-3gr 2 -30zg 2 r 3 (5.22) 

Writing ip' HD {r) — (^^(r) = yields z c = —1/10, gv c = 1/3, and gt c = 1/3, with a 
critical exponent 7 = —1/3, corresponding to the Lee- Yang edge singularity (conformal 
field theory with central charge c(2, 5) = —22/5) coupled to 2D quantum gravity. 

The inclusion of vertices of odd valences does not give any additional multicritical 
singularities. This is why we choose to stick here and in the following to the even case as 
much as possible. 

5.4- Higher genus 

To keep the full fledge of the model, we must keep track of all shifts of indices in ( |5.10| ) . 
This is easily done by still introducing r{z = m/N) = r m , but by also keeping track of 
finite shifts of the index m — > m + a, namely, setting e = 1/N, via r(z + ae) = r m+a . In 
other words, as N — > 00, we still assume that r m becomes a smooth function of z = m/N, 
but keep track of finite index shifts. Solving eq.( |5.1U| ) order by order in 1/N involves 
writing the "genus" expansion 

r( z ) = ^e 2k r^ k \z), (5.23) 

A;>0 

implementing all finite index shifts by the corresponding e shifts of the variable z, and 
solving for the r^'s order by order in e 2 . We finally have to substitute the solution back 
into the free energy ( |5.13|) , with r» = r{i/N). This latter expression must then be expanded 
order by order in e using the Euler-MacLaurin formula. Setting h(x) = (l — x)Log(r(x)/x), 
this gives 

mp- = Jf X> (Jf) = f h ^ dz + f W 1 ) - HO)) + ^(h'(l) - h'(0)) + ... (5.24) 

in which we must also expand r(x) according to ( |5.23| ). The result is the genus expansion 
Fn(V) = ^ N 2 ~ 2h F( h '(V), where F^ h ' is the generating function for graphs of genus h. 
For illustration, in the 4-valent case, we have 

r(z) — gr(z)(r(z + e) + r(z) + r(z — e)) = z (5.25) 
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Writing r(z) = A°\z) + e 2 r^(z) + 0(e 4 ), we find that 

r (1 \z)(l - 6grW(z)) = gr {0 \z)r^" (z) (5.26) 

at order 2 in e, while r^°\z) = (1 — a/1 — 12gz)/(6g), and (1 — Qgr^ ' (z))r(°> (z) = 1, so 
that r^\z) = gr(°\z)r(°> (z)r ( -°> (z). At to order 2 in e, this gives 

i r 1 

F« = — (/V (l) -h' (0))+g dx(l - x)r( 0)/ (x)r (0) "(x) (5.27) 

12 Jo 

where /loO^) = (1 — x)Log(r^ (x) / x) , namely 

FW = ±-Y 9 -3«(4" - ( 2n )) 
24 ^ n v \nj 



n>\ 



9 , 15 2 



(5.28) 



+8 2 



= 8 ^ 




where we have displayed the genus one 4-valent graphs with up to two vertices, together 
with their inverse symmetry factors. 

5.5. Double-scaling limit 

The idea behind the double-scaling limit is to combine the large iV limit and the 
singularity structure of the free energy at all genera into a single scaling function. Let us 
first consider the 4-valent case (|5.25|) . We wish to approach the critical value g = g c = 1/12 



displayed by the planar solution R — (1 — y/1 — 12g)/(6g) at t — 1, at the same time as 
iV — > oo. Setting p(z) = gr(z), we have 

gz = p(z)(l-(p(z + e)+p(z)+p(z-e))) 

(5-29) 
g c = p c {l - 3p c ) 

with g c = 1/12 and p c = 1/6. Subtracting both lines of ( |5.29| ) and expanding up to order 
2 in e yields 

g c -gz = p c (l - 3 Pc ) - p(z)(l - 3p(z)) - e 2 p{z)p"{z) + 0(e 4 ) (5.30) 
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This suggests to introduce rescaled variables and functions g c — gz = a 2 g c y, while p(z) 
p c (l — au(y)), and to expand up to order 2 in a as well: 



a g c y = 3p c a u(y) - e p c a u (y) (5.31) 

where we have noted that dz = —a 2 dy at g = g c . The large iV limit of Sect. 5. 2 is recovered 
by taking e = 0, in which case we are left with u(y) = y/y, another way of expressing the 
planar singularity of the free energy u(y) = y~ 7 , with 7 = —1/2. For non-zero e, all terms 
in ( |5.31| ) will contribute if we take e 2 = a 5 . We then have 



y = u{y) 2 - l -u"(y) (5.32) 

which is nothing but the Painleve I equation. Moreover, the singular part of the free energy 
reads 

F = F sing = N 2 [ dz(l - z)Log(^-)\ sing = N 2 a 5 [ (y - x)u(y)dy (5.33) 

Jo 9 Z Ja- 2 

where g c — g = a 2 g c x. Differentiating twice w.r.t. x yields u(x) = —F"{x). To summarize, 
if we take simultaneously iV — > 00 and g — > g c , but keep the quantity 

Nif^^A=x (5.34) 

fixed, then the singular parts of the free energy at all genera recombine into a single scaling 
function F(x), whose second derivative satisfies the Painleve I differential equation. To 
recover the leading singularity at genus h, we simply have to expand the solution of (|5.32j ) 
at large x as u(x) = J2h>o UhX ^^ 1 ~ 5h * > anc ^ s °l ve the resulting recursion relation for uh- 
This is the so-called double-scaling limit of pure 2D quantum gravity. 

We may repeat this exercise with the multicritical models of Sect. 5. 3, however algebra 
becomes cumbersome. Let us instead look at the scaling limits of the operators P and Q 
acting on the orthogonal polynomials ( |5.5[ ). Let us rescale the orthogonal polynomials p n 
to make them ortho normal, namely set p n = p n /Vhn, so that (|5.7|) (with s n = 0) becomes 
more symmetric 

(Qp)n = Ap n = y/r n+1 p n+1 + y/r^pn-! (5.35) 

or equivalent ly 

Vn,m = \Pmi ^iPnj — y r n+l^m,n+l ~r yfnOm,n—l [p.O'Oj 
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Let us now take the large iV limit. Setting e = 1/N as before, we note that the shift 
operator <S mjn +i, acting on sequences (a m ) can be generated as e ed l dz 5 acting on the 
continuum limit of (a m ), i.e. a function a(z) = a m , for z = m/N. Indeed, one just has to 
write 

^5 m , n+ ia m = a n+1 = e^ a n ~ e e ~* a(z) (5.37) 

m 

Setting again r(z) = r c (l — au(y)), this permits to rewrite Q as 

Q ~ \/r(z) (e e ~^ + e~ 



d 
dz 



dz' 
d 



y/r e (l - au(y))(2 + e 2 -^ + 0(e 4 )) (5.38) 



2v^ - v^(cm - (e-^) 2 + 0(e 4 , a 2 )) 
In the general multicritical case, we must set 

t c -tz = a m+1 t c y (5.39) 

and y = x at z = 1, so that ck ~ a rn+1 dy. The two terms in the r.h.s. of ( |5.38| ) are of 
the same order a provided (ed/dz) 2 = (ea~ m ~ 1 d/dy) 2 is of order a, and we obtain the 
double- scaling condition that e 2 a~ 2rn ~ 2 = a, hence 

N 2 = a" 2 ™ -3 (5.40) 

or equivalently 

N^+3 I _£ j = x (5.41) 

remains fixed while N — > oo and t —>■ t c . Retaining only the coefficient of a, we find that 

A 2 

in the double-scaling limit. This limit is a differential operator, acting on functions of the 
rescaled variable y. 

Let us now turn to P. It will be useful to slightly change the definition of the operator 
P, in the following way 

d\p m (X)e 2 — e 2 p n (A) 

= -y(Pm,V / (Q)&) + (p m ,&) < 5 - 43 ) 

N 



Q^^-2- <y) ( 5 - 42 ) 
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where A is a lower triangular matrix A nm = if n < m. Upon an integration by parts we 

may as well write 

N 

-*n,m J7 - ' \^c )n,m * -™-n,m 

— ~lfV \Q)n,m A nm 



where the matrix A 1 is upper triangular. Eq. ( p.44j ) permits to compute the matrix elements 



of P in terms of those of Q only, by using the first equation when n < m (A nrn = 0) and 
the second one when n > m (A^ m = 0). This can be summarized by the following operator 
relation: 

P = y(nO)+-V"(0)-) (5-45) 

where the index + (resp — ) indicates that we retain only the upper (resp. lower) triangular 
part. In particular, as it is expressed polynomially in terms of Q, P has a finite range, 
namely P n ,m = if \n — m\ > B, B some uniform bound, independent of iV (B depends 
only on the degree of V). This bound ensures that P goes over in the double scaling limit 
to a differential operator of finite degree p, of the form 

1 



aJr c 



: (d p + v 2 dP- z + v 3 d p ~* + ....+ v p ) (5.46) 



to ensure the correct normalization of [P, Q] = 1. From the precise form of P ( 5.45 ), and as 
each derivative d w.r.t. y carries a prefactor ea~ m ~ 1 = \/a, we must have iV x a p < 2 — 1/a, 
which together with the double-scaling condition (|5.40| ) fixes the degree 

deg(P) = 2m + 1 (5.47) 



We must finally write the canonical commutation relation ( |5.6| ) [P, Q] — 1, with the renor- 
malized values P = <i 2m_1 + V2d 2m ~ 2 + .... + vi m and Q = (d 2 — u). Let us introduce the 
square root L of Q, namely the unique pseudo-differential operator 

L = d + J2^d~ l (5-48) 

such that L 2 = Q. This equation is expressed as a triangular system for the £'s, provided 
we normal-order the result by pushing all functions to the left of powers of the differential 
d, by means of the Leibnitz formula d~ l f{y) = J2j>o(~ 1)^ ( l+J ~ )f^(y)d~ l ~^. Let us 
now express P. Solving [P, Q] = rather than [P, Q] = 1 makes no difference as far 
as we only write the equations for the coefficients of positive powers of d: solving these 
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equations precisely allows to express P as a function of Q. As the solution to [P, Q] = 
for a pseudo-differential operator P of degree 2m + 1 is nothing but P = L + , we simply 
have P = (L 2m+1 ) + , where the subscript + indicates that we have retained only the 
differential polynomial part. So far, we have solved all the equations obtained by setting 
to the coefficients of all positive powers of d in [P, Q] = 1. We still have to write the 
d° coefficient. Writing (L 2m+1 )_ = L 2m+1 - (L 2m+l )+ = R m+1 [u]d- x + 0(d~ 2 ), this last 
equation reads simply 

2R m+1 [u]' = l =► 2R m+1 [u] = y (5.49) 

This is nothing but a higher order generalizations of the Painleve I equation, related to 
the so-called KdV hierarchy. From their definition, the "KdV residues" i? m [w] satisfy the 
recursion relation 

R m+1 [u]' = -R m [u]'" ~ \u'Rm[u] - uR m [u]' (5.50) 

obtained by writing (L 2m+1 )_ = ((L 2 ™" 1 ) _Q) _ = (^(L 2 ™- 1 ).)., while the ini- 
tial term reads R\[u] = —u/2. Again, plugging the large y expansion u(y) = 
^2h>o u hy m+1 ~ m+ m ^° eq (|5.49| ) yields a recursion relation for the Uh and gives 



acces to the all genus singular part of the free energy via the relation F s i ng (x)" = —u(x). 
The actual general solution of [P, Q] = 1 involves integration constants which we 
have all set to zero for convenience, hence the most general solution for a degree 2m + 1 
differential operator P reads 

m+l 

2J2nRjM=V (5-51) 

i=i 

for some integration constants \ij. This equation interpolates between the various mat- 
ter critical points [ij = Sj^+x, corresponding to the various multicritical points already 
identified as c(2, 2k + 1) CFT coupled to 2D quantum gravity. From the point of view of 
the m + 1-critical model, the /x's are just dimensionful parameters coupled to the order 
parameters of the theory. 

5.6. Generalization to multi-matrix models 

A large class of multi-matrix models turns out to be solvable by exactly the same 
techniques as those developed in the previous sections for the one- matrix model. It corre- 
sponds to matrices M±, ...,M p with a chain-like interaction, namely involving a quadratic 



form Q a ,b as in (|2.21| ), for which only the elements Q a , a , a = 1,2, ...,p and say Q a ,a+i, 
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a = 1,2, ...,p — 1 are non-vanishing. In this particular case only, the unitary group in- 
tegrations may be disentangled from the eigenvalue integrations for all M's and we may 
still reduce the integral to one over eigenvalues of the different matrices. Once this step is 
performed, the orthogonal polynomial technique is easily adapted and a complete solution 
follows from considering again operators P a and Q a of differentiation w.r.t. or multiplica- 
tion by an eigenvalue of the matrix M a , a = 1, 2, ..., p. Note that the saddle-point technique 
with several matrices is more subtle. 

One is eventually left with solving an equation of the form [Pi, Qi] = 1, the scaling 
function u such that u" = —F being identified with some coefficient of Q\. The remarkable 
fact is that both Pi and Qi remain of uniformly bounded range, the latter depending only 
on the degrees of the potentials for the various matrices. This implies that in a suitable 
double scaling limit where the size of the matrices is sent to infinity and the parameters 
of the potentials go to some (multi-) critical values, the operators Pi and Qi still become 
differential operators of finite degree say p and q, two coprime integers. The resulting 
differential system [Pi,Qi] = 1 governs the all-genus singular part of the free energy of 
the general c(p, q) minimal conformal field theory coupled to 2D quantum gravity. This 
completes the picture of critical behaviors covered by matrix models solvable by orthogonal 
polynomial techniques: it exhausts all minimal CFT's with c < 1, according to the famous 
ADE classification thereof [Qi. 

From a combinatorial point of view, the bijection presented in Sect. 4.1 may be gen- 
eralized to the case of two- matrix models p3fl , and presumably to all cases solvable by 



orthogonal polynomial techniques, which all lead to algebraic systems, henceforth suggest 
tree-like interpretations. 

6. The combinatorics beyond matrix models: geodesic distance in planar 
graphs 

In this section, we return to the bijections between planar graphs and trees to investi- 
gate more refined properties of the discrete random surfaces generated by matrix models, 
involving their intrinsic geometry. In particular, we will derive in a purely combinatorial 
manner sets of closed equations for generating functions of planar graphs with marked 
points at a given geodesic distance, a task still eluding the matrix model description. 



This statement is not completely correct: only the A-type CFT's are covered by the standard 



multi-matrix models. A proposal for D-type CFT's was given in [21], based on D-type generaliza- 



tions of the KP hierarchy [22|, but no direct relation to solvable matrix models was found. Not 



to speak about E-type solutions... 
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6.1. Keeping track of the geodesic distance: the J^-valent case 

Let us return to the bijection between two-leg 4-valent planar graphs and rooted 
blossom trees shown in Sect. 4. Looking at Fig. 6, we see that the bijection allows to keep 
track of the geodesic distance between the two legs, namely the smallest possible number of 
edges of the graph crossed by a curve joining them. Indeed, this distance is nothing but the 
number of edges encompassing the root of the corresponding blossom-tree, when black and 
white leaves are re-connected. Loosely speaking, the geodesic distance between the legs 
corresponds in the blossom tree language to the number of black leaves "in excess" , which 
require encompassing the root to be connected to their white alter ego in counterclockwise 
direction. Let us now derive simple relations for the generating function R n for two-leg 
diagrams with geodesic distance at most n between the legs. To get the more interseting 
generating function G n for graphs with two legs at geodesic distance equal to n, we just 
have to write G n = R n — R n -i- Alternatively, R n can be thought of as the generating 
function for blossom trees with at most n black leaves in excess. As such, it obeys the 
following recursion relation: 



Rr 



1 



mm 




(6.1) 



This is just a refinement of eg. ( |4.1|) in which we have kept track of the maximal numbers 
of excess black leaves. The presence of single black leaves around the vertex connected to 
the root lowers by 1 the maximal number of excess leaves of any object on its left, while 
as each blossom tree has one white leaf in excess, it always absorbs one excess black leaf 
from objects on its left: these two facts are responsible for the shifts of the index n. 

The recursion relation ( |6.1|) holds for all n > provided the term involving R-i is 
dropped. Let us therefore supplement the recursion relation with the initial value R-\ = 0. 
Moreover, the function R n should go over to the function R of (LI) in the limit n — > oo, 
which amounts to suppressing the constraint on the distance between the two legs. R n 
is the unique solution to ( |6.1[ ) such that i?_i = and lim^^oo R n = R. If we are only 
interested in the power series expansion of R n in g, we may solve ( |6.1| ) order by order 
in g, starting with R-\ = at all orders and R n = 1 + O(g) for all n > 0. To any 
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given order in g, the system for the series coefficients is indeed triangular, and moreover 
R — R n = 0(g n+1 ), which guarantees the convergence condition. In the next section, we 
actually display the exact solution R n in a very compact form. 

An important remark is in order. The relation ( |6.1|) is strikingly reminiscent of that 
for the orthogonal polynomials ( |5.11| ) say with g 2 = g§ = 0, g± = g, except that the 
l.h.s. of (|5.11|) is now replaced by 1, and r n = h n /h n -x by R n . One may wonder whether 



eq.( |6.1|) may be derived from some matrix model solution. The answer is not known to 
this day, but the boundary condition that R-± = would mean in matrix model language 
that some norm of orthogonal polynomial must vanish, hence if there is such a matrix 



model formulation, it must be very singular. As to the r.h.s. of (|5.11| ), its similarity with 
that of (|6.1|) suggests to express the rules for the possible subtrees encountered around the 
vertex attached to the root in counterclockwise order in terms of a "Q-operator" acting on 
a formal orthonormal basis \n), (m\n) = b~ m ^ n for m, n > and \n) = for n < 0, via 

Q\n) = \n + 1)+R n \n-1) (6.2) 

The first term is interpreted as the contribution of a single black leaf, while the second 
corresponds to a blossom tree with at most n excess black leaves. Then the r.h.s. of ( |6.1| )is 
nothing but 1 + g(n — l|Q 3 |n). 

6.2. Exact solution 



To solve ( |6.1|) , we use the convergence condition to write R n = R — p n at large n, and 
expand (|6.1|) at first order in p n . This gives the linear recursion relation 

gRipi 1 ^ + pi 1 !,) - p«(l - 4gR) = (6.3) 

This has the characteristic equation 

x + - + 4=— - (6.4) 

x gR 

with R given by (|3.39| ). Picking the solution x with modulus less than 1, we find that 
R n = R(l — Xix n + 0(x 2n )) for some integration constant X\. We may next expand 
R n = R(l — J2j>i AjX jn ), and ( |6.1|) turns into a recursion relation for the coefficients Xj\ 



Ai+i (^' +1 + ^r-^-^)=E X ^+^ ( x% + ^) 



(6.5) 
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solved recursively as 



X±x V 1 — x J 



(1 — x)(l — x l ) J 1 — x 



Picking Ai = x(l — cc)(l — x )A, R n is easily resummed into 

7? _z ? ( 1 -^ +1 )(l-Ax-+ 4 ) 
n (1-Ax™+ 2 )(1-Aa-+ 3 ) { } 

Further imposing the initial condition i?_i = fixes A = 1, hence finally 

(l-x n+1 )(l-x n + 4 ) 
Rn ~ R (l-x^)(l-x^) (6 " 8) 

with |x| < 1 solving (|6.4|) . 

This gives an explicit formula for the generating function of 4-valent two-leg graphs 
with geodesic distance at most n between the legs. In particular, for n = 0, this gives 
the generating function for graphs with the two legs in the same face (also called 1?2 in 
Sect. 3. 3), namely 

R = G = T 2 = R 1 + x2 9 = R ^}~ A9R \ =R-gR 3 (6.9) 

U U 1 + x + x 2 (l-3gR) K ! 



where we have used (|6.4|) and ( |4.1| ) to simplify the result. This is in perfect agreement 
with the matrix model result (|3.41|) . 

6.3. Integrability 

The equation ( |6.1|) is intergable in the classical sense that there exists an "integral 
of motion", namely a conserved quantity f(R n ,R n+ i) =const. which implies flBTl]). More 
precisely, defining 

f(x,y) = xy(l-gx-gy) -x-y (6.10) 

we have 
f(R n , R n +i) — f(Rn-l, Rn) = (R n +i — R n -i) [Rn — 1 — gR n (R n +i + R n + R n -i)) (6.11) 

We deduce that if f(R n , R n +i)is a constant independent of n, then R n obeys (|6.1|) : / 
is an integral of motion of the equation ( |6.1| ) . 

Using /, we may write in a compact way the condition lim n ^ 00 R n = R for solutions 
of ( |6.1| ). Indeed, we simply have to write 

f(Rn, R n+1 ) = f(R, R) = R 2 (1- 2gR) -2R= -{R - gR 3 ) (6.12) 
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All solutions to (|6.12| ) are also solutions of (|6.1|), and they moreover converge to R as 
n — » oo. As an immediate application of (|6.12|) , we may recover Rq, by imposing that 
R-! = 0: 

R = -f(R-u Ro) = R- gR 3 (6.13) 



in agreement with (|6.9| ). 



Rn,A = </> ^-ZAT.Rn (6.14) 



6.4- Fractal dimension 

The advantage of having an exact formula like ( |6.8|) is that we may also extract the 
"fixed area" coefficient R Uj a of g A in R n via the contour integral 

dg 

1mg J 

with R n given by (|6.8|). This gives access to asymptotic properties at large area A. In 
particular, the ratio 

B n = lim ^ (6.15) 

may be taken as a good estimate of the average number of points at a geodesic distance less 
or equal to n from a given point in random 4-valent graphs of infinite area. It is expected 
to behave like 

B n ~ n dp for large n (6.16) 

where dp is the fractal dimension of the random surfaces. Performing in ( |6.14j ) the change 
of variables v = gR, i.e. g = v (1 — 3v), we obtain 

f dv(l-6v) 1 (l-x(v) n+1 )(l-x(v) n+4 ) 

n ' A J 2iir(v(l- 3v)) A + 1 1-3^(1 - x(v) n + 2 )(l - x(v) n + 3 ) [ ' ' 

where we have used R(g(v)) = 1/(1 — 3v) and the expression x = x(v) = (1 — 4v — 

a/1 — 8v + 12f 2 )/(2f). The large A behavior is obtained by a saddle-point approximation, 

as the integral is dominated by the vicinity of v = v c = 1/6, corresponding to the critical 

point g = g c = 1/12, where x — > 1. Making the change of variables v = v c (l + i~rr), 

expanding all terms in powers of 1/y/A and integrating over ^, we finally get the leading 

behavior 

fl2 N ) A (n + Din + 4) 

R nA ~ const. ^4-7 T7 f (140 + 270n + 179n 2 + 50n 3 + 5n 4 ) (6.18) 

A* {n + 2)(n + 3) 

which finally gives the ratio 

„ 3 (n + l)(n + 4) , 9 , A . 3 A , 

B ™ = ^rr^TT^ 140 + 270n + 179n + 50n + 5n ~ ^ n 6 - 19 

280 (n + 2)(n + 3) 56 

hence dj? = 4 is the desired fractal dimension. 
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6.5. Scaling limit: Painleve again! 

A continuum limit may be reached by letting g tend to its critical value g c = 1/12. 
More precisely, we write 

9= T2 (1 ~ e4) => gR=\{l-e 2 ) (6.20) 

from eq.( |3.39| ). In turn, the characteristic equation (fOj) yields 

x = e -^ + 0(e 3 ) a = V6 (6.21) 

As seen from eq.( |6.8|) , a sensible limit is obtained by writing 

n=- (6.22) 

and letting e — > 0. Writing the scaling variable r as r = n/£, we see that e plays the 
role of the inverse of the correlation length £. As we approach the critical point, we have 
£ = e _1 = ((g c — g)/g c ) with a critical exponent v = 1/4, in agreement with v = 1/oIf, 
as expected from general principles. Performing this limit explicitly on the solution ( |6.8| ) 
yields an explicit formula for the continuum partition function T(r) of surfaces with two 
marked points at a geodesic distance larger or equal to r: 

Fir) = lim R ~^ n = -2^Log sinh (y/Z r) = ^ (6.23) 

V ^ e-o e 2 R dr 2 5 VV 2 ! sinh 2 (^f r) V ' 

Upon differentiating w.r.t. r, we obtain the continuum partition function for surfaces with 
two marked points at a geodesic distance equal to r: 

g(r) = -T> (r) = 3V6 ™^f*\ (6.24) 

smh (vfrj 



This reproduces a conjecture |2^] obtained in a transfer matrix formalism of 2D quantum 
gravity. 

Note that the precise form of the scaling function T{r) may alternatively be obtained 
by solving the continuum counterpart of eq.(B7j]). Indeed, writing 



R n = R{l-e 2 F{ne)) (6.25) 

and expanding eq.( |6.1| ) up to order 4 in e, we obtain the following differential equation 

F"(r) - 3F 2 (r) - 6JF( r ) = (6.26) 
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It is easy to check that J-{r) as given by ( |6.23|) is the unique solution of (|6.26| ) with 



boundary conditions T(r) —>■ oo when r — ► and J-"(r) —>■ when r — > oo. Writing 
jF(r) = u(r) — 1, we note that eq.( |6.26| ) turns into 



u 1 - u"/Z = 1 (6.27) 

strikingly reminiscent of the Painleve I equation governing the model's all-genus double- 
scaling limit Q5.32 ), except for the r.h.s. which is now a constant. The function u leading 



to T is simply the unique solution to ( |6.27|) such that u(0 + ) = oo and u(+oo) = 1. 

6. 6. Generalizations 

The results of Sects. 6. 1-6.5 generalize straightforwardly to the case of arbitrary even 
valences. Using again the bijection of Sect. 4. 2, we still have to keep track of excess black 
leaves. Introducing similarly the generating function R n for planar graphs with even 
valences and with two legs at geodesic distance less or equal to n, we get a recursion relation 
by inspecting all configurations of the vertex attached to the root of the corresponding 
blossom trees. We may use the same rules as those found in the 4-valent case ( |6.1|) . 
Going clockwise around the vertex and starting from the root, we may encounter blossom 
trees with up to p excess black leaves or single black leaves. Encountering a black leaf 
decreases the index p of the objects following it clockwise, while encountering a blossom 
subtree increases it by 1. Using the "Q-operator" formalism of Sect. 6.1, namely that 
Q\n) = \n + 1) + R n \n — 1), we get the general recursion relation 

R n = l + J292k(n-l\Q 2k - 1 \n) (6.28) 

fc>i 

to be supplemented with d/2 — 1 initial conditions R-i = R-2 = ...Rd/2-1 — 
(d = deg(V)), and the usual convergence condition lim n ^ 00 R n = R, to the solution 
R of ( f4.2|) . The explicit solution to (|6.28| ) with these boundary conditions was derived in 
||14j| , and involves soliton-like expressions. It allows for investigating the fractal dimension 
for multicritical planar graphs, found to be dp = 2(ra + 1) for the case of Sect. 3. 6 ( |5.21|) , 
and to derive continuum scaling functions for multicritical matter on surfaces with two 
marked points at a fixed geodesic distance r. Writing ( |6.28| ) as 1 = (n — l\V'(Q)\n), we 
use again the trick of adding a weight t per face of the graph, which amounts to replacing 
V(Q) — > V'(tQ), and multiplying by t leaves us with t = (p(t' 2 R n ,t 2 R n ±i, ...). Taking 
the multicritical values for g^k-, and writing t = t c (l — e 2< - m+1 - ) ), we look for solutions of 
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( |6.28|) of the form R n = R(l — ^Tiy = ne)). This gives at order 2(m + 1) in e a differ- 
ential equation for T. Noting that our scaling Ansatz for R n is the same as that for the 
double-scaling limit (r n = r c (l — au(y))) except for the prefactor R = R c (l — e 2 ) we see 
that u(r) = 1 + T{r) satisfies the generalized Painleve equation (5.49), but with a con- 



stant r.h.s. In differentiated form, this corresponds to writing the commutation relation 
[P, Q] = between two differential operators P and Q of the variable r, with respective 
orders 2m + 1 and 2, with Q — d 2 — u. 

The generalization to graphs with arbitrary (even and odd) valences is straightforward, 
as we simply have to use the "Q-operator" formalism in the combinatorial setting. The 
functions S n (resp. R n ) generate planar graphs with one leg (resp. two legs), with the leg 
(resp. second leg) at distance at most n from the external face. The operator Q now acts 
as Q\n) = \n + 1) +S n \n) + Rn\n ~ 1), where the new contribution corresponds to subtrees 
of charge 0, that do not affect the numbers of allowed excess black leaves of their followers. 
We obtain the system of equations 

= (n\V'(Q)\n) 1 = (n - l\V'(Q)\n) (6.29) 

This generalizes presumably to all planar graph enumeration problems for which a 
matrix model treatment is available, using orthogonal polynomials involving a natural Q 
operator, interpreted in the combinatorial setting as describing objects of various charges 
attached to the root vertex of the corresponding blossom trees. We may infer that in the 
general multicritical case of a CFT with central charge c(p, q) < 1, the scaling function for 
surfaces with two marked points at geodesic distance at least r is governed by a differential 
system of the form [P, Q] = 0, P and Q two differential operators of the variable r of 
respective degrees p and q. 

7. Planar graphs as spatial branching processes 

This last section is devoted to a dual approach to that followed so far, in which we 
consider the graphs dual to those contributing to the matrix model free energy, namely 
with prescribed face valences rather than vertex valences. On such a graph, the geodesic 
distance between vertices is the minimal number of edges visited in a path from one to 
the other. We will present bijections between classes of such graphs with a specified origin 
vertex and with a marked vertex at geodesic distance < n, and labeled trees of arbitrary 
valences obeying some specific labeling rules. 

This allows to make the contact with an active field of probability theory dealing with 
spatially branching processes. The following is largely based on refs. |13[] |26| |27[] p8| . 

We first concentrate on the quadrangulations, namely the duals of 4-valent graphs. 
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(a) 



(b) 



Fig. 8: The bijection betwen planar quadrangulations and labeled trees. A 
planar rooted quadrangulation (a) and the natural labeling of its vertices by 
the geodesic distance to the origin vertex of the rooted edge (arrow). The 
confluent faces are shaded. The tree edges are represented in thick black 
lines, and connect all vertices with positive labels. Erasing all but these new 
edges and the vertices they connect leaves us with a labeled tree (b), which 
we root at the vertex corresponding to the end of the rooted edge of the initial 
quadrangulation. Finally, all labels of the tree are shifted by — 1. 

7.1. The dual bijections: labeled trees for planar quadrangulations 

We start with a rooted planar quadrangulation, namely a graph with only 4-valent 
faces (squares), with a marked oriented "root" edge. Let us pick as origin vertex the 
vertex at which the root edge starts. This choice induces a natural labeling of the vertices 
of the graph by their geodesic distance to this origin, itself labeled (see Fig. 8 (a) for an 
example). We then note that only two situations may occur for the labeling of vertices 
around a face, namely 



n+1 



n 



n 



n+1 



n+1 



n+2 



n 



n+1 



(7.1) 
in which cases the faces are respectively called confluent and normal. The confluent faces 
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have been shaded in the example of Fig. 8 (a). We now construct new edges as follows: 




n+1 



n+2 



n 



n+1 



(7.2) 

in each face of the quadr angulation (including the external face, for which the rules are 
reversed). This rule may be summarized by saying that we connect via a new edge all 
the vertices immediately followed clockwise by a vertex with a label one less. These 
edges are readily seen to connect all vertices of the quadrangulation but the origin. Thus, 
erasing all but the new edges and the vertices they connect leaves us with a connected 
labeled tree (see Fig. 8 (b)), which we root at the end vertex of the original rooted edge 
of the quadrangulation, and in which we subtract 1 from all vertices^ In particular, the 
vertex attached to the root has label 0, and all labels are non-negative. Moreover, by the 
construction rules ( [7.2| ), adjacent labels of the tree may differ only by or ±1. Such trees 
are called well-labeled, and are in bijection with the rooted planar quadrandulations. 

The construction rules ( |7.2j ) allow for interpreting the features of the tree in terms of 
the original quadrangulation. Any vertex labeled n — 1 in the tree corresponds to a vertex 
at distance n from the origin in the quadrangulation. From the rules of eq.( [7.2|) , we see 
that any marked edge n —>■ n + 1 of the quadrangulation corresponds marking an edge of 
the tree adjacent to a vertex labeled n. This in turn may be viewed as the rooting of the 
tree at a vertex labeled n (the above bijection uses this fact for n = 0). 

We next define rooted well-labeled trees as rooted labeled trees, with non-negative 
integer vertex labels, and such that the root vertex has label n. Let R n be the generating 
function for such objects, with a weight g per edge. According to the above bijection, 
the generating function for rooted planar quadrangulations with a weight g per face is 
simply Rq. If, instead of rooting the well-labeled tree at the end vertex of the initial 
quadrangulation, we had chosen to root it elsewhere, typically at another vertex of the 
tree say labeled n, the resulting rooted well-labeled tree would satisfy the extra condition 
that the label occurs at least once in the tree. The generating function for such an 



This is just a technical trick to make the precise contact with the generating function R n of 
Sects. 4 and 6. The reader will have to remember to add up one to each vertex label of the tree to 
recover its geodesic distance from the origin in the quadrangulation. 
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object is nothing but G n = R n — R n -i- In terms of the original quadrangulation, this is 
nothing but the generating function of quadrangulations with an origin vertex and with 
a marked edge n —>■ n + 1 w.r.t. this origin. So R n is the generating function for planar 
quadrangulations with an origin and with a marked edge m—>-m + l,m<n, and a weight 
g per face. 

The definition of R n allows to derive a recursion relation of the form 

1 ~~ g(Rn+l + Rn + Rn-l) 

where we simply express the labeling rule that the root vertex labeled n may be adjacent 
to any number of vertices labeled n, n + 1 or n — 1, themselves roots of other well-labeled 
trees. Moreover, for ( |7.3| ) to also make sense at n = we must set R-i = 0. Removing 
the constraint that m < n by sending n — > oo leaves us with the generating function R 
for quadrangulations with an origin and a marked edge, wich also generates the rooted 
quadrangulations with a marked vertex, and should satisfy the relation 

with R = 1 + O(g). We conclude that the functions R and R n coincide with those intro- 
duced in Sects. 4.1 and 6.1. 

So we have found another (dual) combinatorial interpretation for the exact solutions 



7.2. Application I: average numbers of edges and vertices at distance n from a vertex in 
quadrangulations 

A direct application of this new interpretation of R n concerns properties of large 
random quadrangulations viewed from their origin. For instance, the average (e n )^ of the 
number of edges n — > n + 1 in a quadrangulation with an origin and with say A faces is 
given by 

\ e n/ Rn,A ~ Rn — 1,A /-, r\ 

(e ) Rq,a 

with R n) A as in ( |6.14| ). Again, this is readily computed in the limit A — > oo, where we first 
note (eo) — > 4 by Euler's relation, and then use a saddle point method just like in ( |6.18| ), 
resulting in 

6 (n 2 + An + 2)(5n 4 + 40n 3 + 117n 2 + 148n + 70) 
(Cn) "35 (n + l)(n + 2)(n + 3) ( } 
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This goes as 6n 3 /7 for large n, which confirms the value du = 4 for the fractal dimension, 
as (e n ) sin d/dnn dF ~ n dp ~ 1 . 

We may also obtain the average number of vertices at geodesic distance n from the 
origin, by noting that the corresponding generating function is that of unrooted well- 
labeled trees with at least a label and a marked vertex with label n — 1. Abandoning 
the condition that a label should occur, and decomposing the tree according to the 
environment of the marked vertex with label n — 1 results in the generating function 

oo u 

K n - t = J2 ^j( R n + Rn-i + R n -2) k = -Log(l - g(R n + R n - X + R n -2)) = Lo#(i? n _i) 

k=i 

(7.7) 

where we have incorporated the symmetry factor 1/k when the vertex has valence k. 

Finally, the generating function for quadr angulations with an origin and a marked vertex 

at distance n is 

V n = K n _ t - K n _ 2 = Log (jj^) (7-8) 

for n > 2 and Logi?o for n = 1, while of course Vq = 1. Therefore the average number of 
vertices at distance n from the origin in a quadrangulation of area A is given by 

Ma = Log (^=^) (7.9) 

easily derived in the large A limit: 

(v n ) = -^(( n + !)( 5 ^ 2 + 10n + 2 ) + <W) (7-10) 

This goes as 3n 3 /7 for large n, also in agreement with dp = 4. 

Note that eqs.([7.7|) ( |778|) also allow to interpret Logi? n _i as the generating function 
for quadrangulations with an origin and a marked vertex at distance m < n. In the limit 
n — » oo, the function LogR therefore generates the quadrangulations with two marked 
vertices. In the dual formulation, this corresponds to 4-valent planar graphs with two 
marked faces: this gives a purely combinatorial derivation in the 4-valent case of the 
formula ( |5.19| ) obtained above in the matrix model language. 
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7.3. Application II: local environment of a vertex in quadrangulations 

Another application of this new graph interpretation of R n concerns the local envi- 
ronment of the origin. Assume we wish to keep track of the numbers of vertices at some 
finite distances p + 1 from the origin, and edges labeled q — » q + 1 for some specific p'sand 
g's, both less or equal to some given k. Then a way to do it is to add extra weights, say 
p p per vertex labeled p in the corresponding well-labeled tree and a p per edge adjacent to 
a vertex labeled p of the well-labeled tree. Indeed, as explained in the previous section, 
this amounts to adding a weight p p per vertex labeled p + 1 in the quadrangulation, and a 
weight o p per edge p — > p + 1 in the quadrangulation. This turns the equation ( ff73|) into 
a new set of equations 

Rn = , ^7 5—^- ^— 5 r, n = 0,l,2,...,fc+l (7.11) 



with pfc+i = CTfc+i = 1, while i? n satisfies ( |7.3| ) for all n > k + 2. This is slightly simplified 
by introducing Z n = a n R n (with a n = 1 for n > k + 1), as we are left with 



Zljl T T" , Tl U, 1 , Jin . . . , K 

l — go~n{Z n +i + Z n + Z n -i) 
Z n = z t^ r, n = k + l,k + 2,... 



(7.12) 



Solving such a system seems quite difficult in general, but we may use the integral of 
motion (|6.10| ) to replace the infinite set of equations on the second line of ( [7. 12|) (and the 



convergence condition of Z n to R) , by simply the conserved quantity 

f(Z k ,Z k+l ) = f(R,R) (7.13) 

Together with the first line of ( |7. 12|) , this gives a system of k + 2 algebraic relations for 
the functions Zq, Z\, ..., Z^+i, which completely determines them order by order in g. As 
an example, let us compute in the case k = the generating function including a weight 
Po — P P er vertex labeled in the trees and o~q = a per edge incident to a vertex labeled 
in the trees. (This in turn corresponds in the quadrangulations to a weight p per vertex 
labeled 1, i.e. per nearest neighbor of the origin, and a weight a per edge — > 1.) We get 
the system: 

Z ° = 1 Z r^V Z » Z ^ 1 ~ 9( Z o + Z i) -Z -Z 1 = f(R, R) = gR 3 -R (7.14) 

1 - ga{Z + Z x ) 
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which upon eliminating Z\ and reinstating R = Zq/ct, boils down to 

(Ro -p)(l + R - ga 2 Rl -p)- aR (Ro -p + gR(l - gR 2 )) + g<? 3 R 3 = (7.15) 

for the generating function Rq for rooted quadrangulations with weights p per neighboring 
vertex of the origin and a per edge adjacent to the origin. Rq = Ro(g\p, a) is the unique 
solution to (J7.15Q such that Rq = p + O(g). Note that we recover R = R — gR 3 of 



( |6.13 ) when p = a = 1. As the rooting of the quadrangulation is itself a choice of an edge 
adjacent to the origin, we may express the corresponding generating function for "unrooted 
" quadrangulations, namely with just an origin vertex, as 

r ds 
T (g\p 1 a)= / —Ro(g\p,s) (7.16) 

Jo s 

simply expressing the rooting of the quadrangulation as ad a T Q = Rq. The statistical 
average over quadrangulations of area A of p Nl a No1 (N\ the number of neighboring vertices 
of the origin, Nqi the number edges adjacent to the origin) finally reads 



r ,^(i,i) fi%Ro tA (i, a ) 



{p N la N 0l)A = ^,A^uj = Jo , ^,w, (717) 



where as usual Vq^P, <t) (resp. Rq^Pi s)) denotes the coefficient of g A in r (<?|p, cr) (resp. 
Ro(g\p, s)). The limit \im.A^oo(p Nl o~ Noi )A = T may again be extracted by a saddle-point 
expansion. After some algebra, we find 

er(r + i)(r + 3) - er(2r(i + 4r + r 2 ) + 3p(r + i) 2 (r + 2)) = o (7.18) 

and r is uniquely determined by the condition T = 1 for a = p = 1. For instance, when 
a = 1, we get 

in which we read the probability P(n) = (3/16) n ( ™) for a vertex to have n neighboring 
vertices in an infinite quadrangulation. Similarly, taking p = 1, we get 



which generates the probabilities to have n edges adjacent to a vertex in an infinite quad- 
rangulation. We may also derive the generating function for the conditional probabilities of 
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having n nearest neighboring vertices, given that there is no multiple edge connecting them 
to the origin, by simply taking T(p = t/cr, a) and letting a — ► 0, which indeed suppresses 
all contributions from multiply connected vertices. This gives 



n «=i™ r (H = \& 2 (7 - 21) 

For instance, the probability that a given vertex have no multiple neighbors in an infinite 
quadrangulation is 

n(l) = v / 7-2 (7.22) 

7.4- Spatial branching processes 

We have seen so far how the information on the geodesic distance from the origin in a 
rooted planar quadrangulation may be coded by rooted well-labeled trees. The latter give 
rise to natural examples of so-called spatially branching processes, in the context of which 
quantities like R n correspond to certain probabilities. 

A spatial branching process consists of two data. First we have a monoparental 
population, whose genealogy is described by a rooted tree, the root corresponding to the 
common ancestor. A standard measure on these trees attaches the probability (1 — p)p k 
for any vertex to have k descendents. The second data is a labeling of the vertices of the 
tree by positions say on the integer line n£l Here, we add the rule of the "possessive 
ancestor" that his children must be at close enough positions from his (namely differing 
by or ±1). Let E(T) denote the probability of extinction of the population at generation 
T, then we have the recursion relation 

En{T) = 1 - § (E n+1 (T - 1) + E n (T - 1) + E n _ x (T - 1)) (7 ' 23) 

Letting T — > oo, we see that the extinction probability E n = lim.T—Kx> E n (T) obeys 



the same equation as R n (|7.3| ) upon some rescaling, and we find that E n = (1 — 
p)R n ( g = 7 ] , in the case of positions restricted to lie in a half-line (with a "wall" at 
the origin). Without this restriction, the problem becomes translationally invariant and 
E n = E = (1 — p)R (g = p ^ ~ p ' j . Note that the critical point g = g c = 1/12 corresponds 
here to the critical probability p = p c = 1/2. 

In this new setting, we may ask different questions, such as what is the probability 
for the process to escape from a given interval, say [0, L] . Once translated back into R n 
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terms, this amounts to still imposing the recursion relation (|7.3| ), but changing boundary 
conditions into 

#-1=0 and R L+1 = (7.24) 

The escape probability from the interval reads then 

S n = l- (l-p)R n (g = P ^~ P) ) = (l-p)(R-R n ) (7.25) 



The equation ( |7.3| ) with the boundary conditions ( [7.24| ) still admits an exact solution 
expressed by means of the Jacobi #i function 

01 (z) = 2zsin(7rz) Y[(l - 2q J cos(2ttz) + q 2j ) (7.26) 

The solution reads R n — R n , with 



R{ L) = R UnU n+ 3 
W n +lW n +2 

n + l 



u n = $x 



L + 5 



(7.27) 



guaranteeing that the boundary conditions ( |7.24j) are satisfied, and where the nome q still 
has to be fixed. The main recursion relation (|7.3| ) reduces to a quartic equation for the 

Mn's: 

u n u n+1 u n+2 u n+3 = —u 2 n+1 u 2 n+2 + gR(u n -iul +2 u n+3 + u 2 n u 2 n+3 + u n u 2 n+1 u n+ ±) (7.28) 



and the latter is satisfied by ( |7.27|) provided we take 



_ B x (a) 6i (2a) f9[ (a) 16[(2a 

R = 4 



9 



d[(0)9i(3a) V^i(a) 2^(2a) 

0' 1 (O) 2 1 (3a) (7-29) 



18*(a)%(2a)(g$-ig^ 



for a = 1/(L + 5). The identity ( |7.28| ) is proved typically by showing that both sides have 
the same transformations under n — > n + L + 5 and n —* n + (L + 5)/(2i7r)Logg, and that 
moreover they have the same zeros, this latter condition amounting to ( |7.29| ). 

The elliptic solution R n may be interpreted terms of bounded graphs as follows. The 
quantity G n = Rn — R n -i is the generating function for quadrangulations with an 
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origin and a marked edge n — > n + 1, which are moreover bounded in the sense that all 
vertices are distant by at most L + 1 from the origin. 

Taking again the continuum scaling limit of the model leads to the probabilists' Inte- 



grated SuperBrownian Excursions (ISE), here in one dimension [29]. The scaling function U 
obtained from R n = R c {l — e 2 U) in the limit fl6.20p , while moreover r = ne and A = (L + 5)e 
are kept fixed, reads: 

U{r) =2p(z\u,uj') (7.30) 

where p is the Weierstrass function (p = — S^Log^i), with half-periods u> = A/2 and u/, 
related via the condition that the second invariant g2(w,uj') = 3. 

7. 5. Generalizations 

We have so far only discussed quadr angulations and their relations to spatial branching 
processes (see also [[JO] ||31[]). All of the above generalizes to rooted planar graphs with 
arbitrary even face valences. These are in bijection with rooted well-labeled trees with more 
involved labeling rules, also called well-labeled mobiles |2£|. This allows for a generalization 



of spatial branching processes, possessing these labeling rules. As we already know that 
these objects have an interesting variety of multicritical behaviors, this should turn into 
multicritical generalizations of the ISE. 



In [28], the general case covered by two-matrix models is treated as well, and seen to 
generate Eulerian (i.e. vertex-bicolored) planar graphs. The latter contain as a particular 
case the gravitational Ising model, and in principle allow for reaching any c(p, q) CFT 
coupled to 2D quantum gravity. These will lead presumably to interesting generalizations 
of the ISE. 



8. Conclusion 

In these lectures we have tried to cover various aspects of discrete 2D quantum gravity, 
namely of statistical matter models defined on random graphs of given topology. 

The matrix model approach, when solvable, gives exact recursion relations between 
quantities eventually leading to compact expressions for the genus expansion of the free 
energy of the models. We have further investigated the so-called double scaling limit in 
which both matter and space degrees of freedom become critical, allowing for instance to 
define and compute a scaling function summarizing the leading singularities of the free 
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energy at all genera, as a function of the renormalized cosmological constant x. The final 
general result takes the form 

[P,Q] = 1, P = d p +v 2 d p - 2 + ... + v pi Q = d q + u 2 d q - 2 + ... + u q (8.1) 

with d = d/dx, all v's and u's functions of x, and u 2 proportional to F" , the second 
derivative of the singular part of the all-genus free energy w.r.t x. 

The combinatorial approach, when bijections with trees are available, also gives exact 
recursion relations between basic generating functions which can be interpreted in terms 
of planar graph counting, while keeping track of the geodesic distance between marked 
points. The expressions for the solutions are completely explicit, allowing for taking a 
scaling limit, describing the free energy for random surfaces with marked points at a 
renormalized geodesic distance r. We may write the general result for this scaling free 
energy in the form 

[P,Q] = 0, P = d p + v 2 d p ~ 2 + ... + v p , Q = d q + u 2 d q - 2 + ... + u q (8.2) 

with d = d/dr, all it's and v's functions of r, and u 2 proportional to the scaling two-point 
function for surfaces with two marked points at geodesic distance > r. 

Remarkably, in all cases solved so far, the exactly solvable geodesic distance problems 
for planar graphs all correspond to cases where a matrix model solvable by orthogonal 
polynomials is available. It seems therefore that the bijections with trees exactly parallel 
the orthogonal polynomial solutions. More precisely, we have observed that a similar ab- 
stract "Q-operator" could be introduced in both cases, one of them describing the possible 
subtrees one can encounter when going counterclockwise around a vertex of a blossom 
tree, the other describing the multiplication by an eigenvalue A on the basis of orthogonal 
polynomials. 

The two apparently unrelated results (|8.1|) and (|8.2|) show that something deeper 



happens here, that deserves to be better understood. One may imagine that there must 
exist a more general structure which would unify and combine the notions of genus and 
geodesic distance, and give for instance closed equations for scaling functions of both x 
and r. To reach this, one should first be able to control geodesic distances in higher genus 
as well, by generalizing the tree bijection techniques explained here only in the planar 
case. Another possibility could be that matrix models as we know them today may still 
be only part of a more general setting. Some generalizations of matrix models involving 
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integration of eigenvalues over contours (or linear combinations thereof) in the complex 
plane may be the correct answer, and relate to the intrinsic geometry of graphs once 
interpreted combinatorially. 

Finally, it is interesting to notice that no continuum field theoretical representation of 
geodesic distance dependence of random surfaces has been found yet, although 2D quantum 
gravity is now well understood in terms of the coupling of CFT to the Liouville field theory 
. The simplicity of the results found here for the various scaling functions comes as a 
surprise in that respect. Field theory probably still has some way to go before explaining 
the purely combinatorial results shown here. 
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